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ABSTRACT 

A  Study  of  Magnetic  Activity  Effects  on  the  Thermospheric 
Winds  in  the  Low  Latitude  Ionosphere 

by 

Ricardo  Cruz  Davila,  Master  of  Science 
Utah  State  University,  1994 

Major  Professor:  Dr.  Bela  G.  Fejer 
Department:  Physics 

The  purpose  of  this  thesis  is  to  examine  the  effects  of  magnetic  activity  on  the 
low  latitude  F- region  thermospheric  winds.  The  F-region  (120-1600  km)  is  a  partially 
ionized  medium  where  O*  and  O  are  the  major  ion  and  neutral  species,  respectively. 
The  thermospheric  winds  at  these  altitudes  are  driven  primarily  by  pressure  gradient 
forces  resulting  from  the  solar  heating  during  the  day  and  cooling  at  night.  For  this 
study,  we  use  measured  Fabry-Perot  Interferometer  (FPI)  winds  at  Arequipa  (16.5°S, 
71.5°W)  and  measured  FPI  and  Incoherent  Scatter  Radar  (ISR)  winds  at  Arecibo 
(18.6°N,  66.8°W). 

Previous  wind  studies  at  Arequipa  and  Arecibo  concentrated  on  the 
climatological  wind  patterns  highlighting  solar  cycle  effects  and  seasonal  variations; 
however,  these  studies  did  not  address  the  effects  of  magnetically  disturbed  conditions 
on  the  seasonal  averaged  winds.  To  properly  investigate  storm  time  effects  on  the 
neutral  winds,  we  must  first  investigate  solar  cycle  effects  on  the  seasonal  averages 
during  magnetically  quiet  (Kp  <  3)  conditions.  This  study  will  include  a  detailed  analysis 
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of  solar  cycle  effects  on  the  seasonally  averaged  winds  for  Arequipa  and  Arecibo.  In 
addition  to  the  wind  averages,  we  used  cubic  splines  to  fit  the  average  wind  profiles  and 
to  provide  better  comparisons  with  modeled  results.  We  also  performed  a  study  on  the 
airglow  emission  heights  using  both  Jicamarca  and  Arecibo  electron  density  profiles. 
This  established  the  height  which  we  will  use  to  compare  our  experimental  data  with  the 
model  winds. 

To  investigate  magnetic  activity  effects  on  the  FPI  and  ISR  winds,  we  used  three 
magnetic  activity  cases  which  cover  all  storm  time  scenarios.  These  magnetic  activity 
cases  are  the  extended  quiet,  short-term  disturbed,  and  extended  disturbed  conditions. 
The  first  case,  the  extended  quiet,  is  the  condition  where  the  previous  and  short  term 
magnetic  activity  is  quiet  (12  hour  Kp  <  3  and  the  Kp  <  3).  The  short-term  disturbed 
case  is  defined  for  the  condition  where  the  previous  magnetic  activity  is  quiet  (12  hour 
Kp  <  3)  and  then  becomes  disturbed  (Kp  >  3).  Last,  we  considered  the  case  where 
previous  and  short-term  magnetic  activity  are  disturbed  (12  hour  Kp  >  3  and  the  Kp  >  3). 

Our  last  objective  is  to  use  our  data  to  validate  the  predictions  from  the 
Thermosphere/Ionosphere  General  Circulation  model  (TIGCM93)  and  the  Horizontal 
Wind  Model  (HWM93).  This  study  should  further  our  understanding  of  the  physical 
processes  which  produce  the  low  latitude  quiet  and  disturbed  winds.  The  TIGCM93  is 
a  first  principal  model  and  the  HWM93  is  an  empirical  model  based  on  ground-based 
and  satellite  measurements.  The  main  advantage  of  using  the  TIGCM93  is  the  ease  of 
studying  the  dynamics  of  ionospheric  phenomena  by  simply  changing  various  model 
inputs,  while  the  HWM93  allows  us  a  comparison  between  our  experimental  wind  data 
sets  with  the  established  climatology  of  the  winds  over  Arequipa  and  Arecibo. 
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CHAPTER  1 


INTRODUCTION 

The  ionosphere  is  defined  as  the  tenuous  gas  of  electrons  and  ions  that  envelop 
the  earth  from  60  km  to  about  600  -  800  km.  Early  radio  experiments  conducted  in  the 
1930s  indicated  the  ionosphere  could  be  subdivided  into  regions  based  on  the  electron 
and  ion  densities.  Today,  the  standard  designations  for  these  regions  of  the  ionosphere 
are  the  D,  E,  Fi,  F2,  and  topside  ionosphere.  Typically  the  Fr  and  F2-regions  are 
referred  to  as  simply  the  F-region.  The  D-region  (60  -  90  km)  is  present  only  during  the 
day  and  is  the  one  region  that  contains  negative  ions.  The  E-region  (90  -  140  km)  and 
Fi -region  (140  -  200  km)  consist  of  molecular  ions,  while  the  F2-region  (200  -  400  km) 
consists  of  atomic  ions  with  being  the  dominant  ion  species.  In  addition  to  ion 
composition,  the  E-  and  F-regions  exhibit  significant  electrodynamic  effects.  Tidal 
oscillations  generate  electric  fields  and  currents  in  the  E-region  while  neutral  wind 
dynamics  drive  the  electrodynamics  in  the  F-region.  The  topside  ionosphere  is  the  name 
given  to  the  ionosphere  above  ~  400  km.  This  region  is  generally  characterized  by  an 
exponential  decrease  in  the  electron  density  and  a  gradual  transition  from  the  ions  to 
the  lU  ions  above  ~  800  km.  Figure  1  shows  the  typical  altitude  profile  of  the 
ionospheric  regions  with  respect  to  the  electron  density.  The  approximate  breakdown 
of  these  regions  by  altitude,  daytime  electron  density,  and  ion  composition  is  shown  in 
TABLE  1.  A  more  detailed  description  of  the  ionosphere’s  regions  and  dynamics  can 
be  found  in  a  number  of  textbooks  [e.g.,  Chamberlain,  1987;  Hargreaves,  1992; 
Kelley,  1989;  Tohmatsu,  1990]. 


2 


Figure  1.  Typical  middle  latitude  daytime  (solid  line)  and  nighttime  (dashed  line) 
electron  density  profiles  and  ionospheric  regions  up  to  600  km  for  moderate  solar  flux 
(^>  =  150)  conditions  and  quiet  geomagnetic  (Ap  =  4)  conditions.  The  electron  density 
profiles  were  generated  using  the  International  Reference  Ionosphere  model  [Bilitza, 
1989]. 
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TABLE  1.  Ionospheric  Regions  by  Altitude,  Major  Ion,  and  Day/Night  Electron  Density. 


Region 

Altitude  Range 
(km) 

Major  Ion 
Components 

Typical  Electron  Densities 
(cm'^) 

D 

60-90 

No^  o; 

1.5  X  lO'^  (day) 
absent  at  night 

E 

90  - 140 

b 

+  cs 

O 

1.5  X  10^  (day) 

<  1.0  X  10"^  (night) 

Fi 

140  -  200 

0^  NO^ 

2.5  X  10^  (day) 
absent  at  night 

F2 

200  -  400 

O^N^ 

1.0  X  10^  (day) 

1.0  X  10^  (night) 

Besides  the  electron  density  profile,  the  atmosphere  can  also  be  classified  by  the 
temperature  gradient  (see  Figure  2).  The  troposphere  (0-12  km)  is  the  region  where 
the  temperature  gradient  decreases  by  approximately  7  K/km.  The  boundary  between 
the  troposphere  and  the  next  region,  the  stratosphere,  is  called  the  tropopause  (~  12 
km).  In  the  stratosphere  (12  -  55  km),  the  temperature  increases  with  height  because 
ozone  absorbs  virtually  all  the  incoming  solar  extreme  ultraviolet  (EUV)  radiation.  The 
stratopause  separates  the  stratosphere  from  the  next  region,  the  mesosphere.  The 
mesosphere  (55  km  -  90  km)  is  characterized  by  a  decreasing  temperature  with  height 
and  contains  the  coldest  temperature  (~  180  K)  at  its  upper  boundary,  the  mesopause. 
The  mesopause  is  the  lower  boundary  of  the  final  temperature  region,  the  thermosphere 
(>  90  km).  The  thermosphere’s  temperature  increases  with  height  since  atomic 
oxygen,  O,  absorbs  solar  EUV  radiation.  The  temperatures  in  the  thermosphere  can 
reach  as  high  as  1300  K  above  300  km. 
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Figure  2.  Typical  temperature  gradient  profile  and  heights  for  the  atmospheric 
temperature  regimes,  liie  temperature  profiles,  provided  by  the  MSISE90  model 
[Hedin,  1991],  for  low  solar  flux  (O  =  75)  and  high  solar  flux  (O  =  200)  conditions. 
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There  are  significant  latitudinal  effects  on  the  ionosphere  density  distribution  and 
dynamics.  At  high  latitudes  the  magnetosphere-ionosphere  coupling  is  the  dominant 
force  driving  the  neutral  and  electromagnetic  dynamics.  The  precipitation  of  particles  in 
the  auroral  oval  directly  affects  the  electron  density  distribution  and  can  be  the  major 
source  of  ionization.  At  middle  latitudes,  the  height  of  the  F-region  electron  density 
peak  is  primarily  affected  by  the  neutral  wind,  diffusion,  recombination,  and  downward 
flow  of  plasma  from  the  ionospheric  heights  above  1200  km.  The  middle  latitude 
maximum  electron  density  concentration  occurs  at  the  height  where  the  downward 
diffusion  and  electron  loss  by  recombination  are  comparable.  The  primary  source  of 
ionization  in  this  region  is  by  the  solar  EUV.  At  low  latitudes,  the  magnetic  field  is 
nearly  horizontal,  which  results  in  unique  transport  effects  associated  with  vertical 
ExB  drifts.  In  addition  to  the  vertical  drifts,  the  electron  density  is  influenced  by 
interhemispehric  flows  caused  by  the  neutral  winds. 

The  focus  of  this  thesis  will  be  the  low  latitude  F-region,  a  complex  and  dynamic 
region  of  the  ionosphere  where  a  variety  of  processes  affect  the  plasma  concentration. 
In  this  region,  the  predominant  neutral  constituent  is  atomic  oxygen,  which  is  ionized  by 
solar  radiation  of  wavelengths  less  than  911  A.  However,  as  altitude  increases  in  the  F- 
region,  the  plasma  density  distribution  is  determined  more  by  transport  processes  such 
as  ambipolar  diffusion,  neutral  wind  drag,  and  electrodynamics  rather  than  local 
ionization.  Plasma  is  transported  along  the  magnetic  field  lines  by  diffusion  and  the 
neutral  wind  effects  and  across  the  magnetic  field  by  electrodynamic  forces. 


6 


For  this  work,  we  study  the  dynamics  of  the  low  latitude  neutral  winds.  Below 
120  km,  the  lower  boundary  of  the  thermosphere,  the  zonal  winds  are  weakest  at  the 
tropopause  (~  12  km)  and  the  mesopause  (~  85  km).  The  strongest  winds  below  120 
km  occurs  near  the  stratopause  (~  55  km).  The  winds  in  the  stratosphere  and 
mesosphere  result  from  the  balance  between  the  pressure  gradient  and  Coriolis  forces. 
Above  120  km,  the  wind  speeds  gradually  increase  and  the  pressure  gradient  force 
becomes  the  primary  forcing  term.  The  pressure  gradient  force  results  from  the  in  situ 
heating  of  the  upper  atmosphere  by  the  absorption  of  EUV  by  atomic  oxygen.  Besides 
pressure  gradients,  the  thermospheric  winds  can  be  influenced  by  geomagnetic 
disturbances  generated  at  the  high  latitudes. 

The  effects  of  geomagnetic  disturbances  (or  storm  time  effects)  on  the 
thermospheric  winds  have  been  observed  at  all  latitudes  and  modeled  using 
sophisticated  theoretical  circulation  models  [e.g.,  Richmond  and  Matsushita,  1975; 
Richmond,  1979;  Blanc  and  Richmond,  1980;  Burnside  et  al,  1991;  Forbes  et  al., 
1993;  Buonsanto  and  Foster,  1993].  Several  cases  of  high  latitude  disturbance  effects 
on  the  thermospheric  dynamics  at  low  latitudes  have  also  been  documented  [e.g., 
Burnside  et  al,  1991;  Fejer,  1993;  Buonsanto  and  Foster,  1993].  One  such  physical 
disturbance  described  in  the  literature  is  the  generation  of  large-amplitude  gravity  waves 
that  propagate  towards  low  latitudes  during  storm  time  conditions.  These  waves  have 
been  shown  to  be  instrumental  in  affecting  both  the  neutral  atmosphere  and  low  latitude 
electrodynamics  [Burnside  et  al,  1991;  Fejer,  1993]. 
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Global  numerical  models  are  a  necessary  and  valuable  tool  for  enhancing  our 
understanding  of  storm  time  processes.  We  use  the  National  Center  for  Atmospheric 
Research  (NCAR)  thermosphere-ionosphere  general  circulation  model  (TIGCM93) 
[Roble  et  ah,  1988;  Fesen  et  al,  1993]  and  the  empirical  horizontal  wind  model 
(HWM93)  [Hedin  et  al.,  1991]  to  assess  the  current  abihty  to  model  these  high  latitude 
disturbances  on  the  low  latitude  neutral  winds.  The  NCAR  TIGCM93  solves  the  full 
nonlinear  primitive  equations  for  eastward  and  northward  momenmm,  thermodynamics, 
continuity,  and  hydrostatics  for  24  constant  pressure  levels  from  97  to  500  km  [Fesen  et 
al,  1993]  while  the  HWM93  uses  wind  data  measured  by  satellite  and  ground-based 
incoherent  scatter  radar  and  Fabry-Perot  optical  interferometers.  The  TIGCM93 
provides  insight  into  what  factors  influence  low  latitude  ionosphere  for  quiet  and 
disturbed  conditions  while  the  HWM93  describes  the  average  wind  behavior  for  a 
selected  altitude,  solar  flux  ,  and  geomagnetic  condition. 

The  primary  focus  of  this  thesis  research  is  to  examine  the  F-region  neutral  wind 
response  to  storm  time  events.  The  data  used  for  this  study  are  winds  measured  by 
Fabry-Perot  Interferometers  (FPI)  at  Arecibo  and  Arequipa,  derived  neutral  winds  from 
Incoherent  Scatter  Radar  (ISR)  at  Arecibo,  and  electron  density  profiles  collected  at 
Jicamarca,  Peru  and  Arecibo,  Puerto  Rico.  Solar  activity  is  determined  by  the  10.7-cm 
solar  flux  index  ($)  for  the  day  of  the  FPI  neutral  wind  measurement,  ISR  derived 
neutral  wind,  and  measured  electron  density.  Magnetic  activity  is  determined  by  the  Kp 
index  for  the  given  observation.  Quiet  and  dismrbed  magnetic  activity  conditions  are 
given  by  Kp  <  3  and  Kp  >  3,  respectively. 
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Chapters  2  and  3  provide  a  brief  introduction  to  the  morphology  of  the  low 
latitude  ionosphere,  upper  atmospheric  wind  system,  and  current  atmospheric  models 
for  investigating  the  wind  climatology  and  physics.  Chapter  4  describes  in  detail  the 
wind  measuring  techniques  and  the  data  analysis  used  for  averaging  the  quiet  and 
disturbed  neutral  winds.  Chapters  5  and  6  provide  the  results  of  the  response  of  the 
neutral  wind  to  quiet  and  disturbed  conditions,  respectively.  These  chapters  also 
provide  comparisons  of  the  measured  winds  to  the  TIGCM93  and  HWM93  predicted 
results.  These  comparisons  help  determine  the  ability  to  model  the  low  latitude  wind 
response  to  quiet  and  disturbed  conditions.  Finally,  chapter  7  contains  the  conclusions 
and  suggestions  for  future  neutral  wind  studies. 
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CHAPTER  2 

THE  LOW  LATITUDE  IONOSPHERE  THERMOSPHERIC  WINDS 

2.1.  INTRODUCTION 

The  focus  of  this  thesis  is  the  low  latitude  thermospheric  neutral  winds.  Before 
investigating  the  neutral  wind  response  to  magnetically  quiet  and  disturbed  geomagnetic 
conditions,  we  must  first  establish  the  morphology  of  the  low  latitude  ionosphere.  This 
chapter  describes  the  unique  characteristics  of  the  low  latitude  ionosphere,  provides  the 
physical  description  of  the  upper  atmospheric  winds,  and  sununarizes  the  low  latitude 
thermospheric  wind  studies. 

2.2.  LOW  LATITUDE  IONOSPHERE 

For  this  study,  we  use  measured  and  modeled  data  from  the  American  sector  (~ 
70"W  geographic)  low  latitude  ionosphere.  The  American  sector  low  latitude 
ionosphere  includes  the  ionospheric  regions  over  Arequipa,  Peru  (16.5°S,  7L5°W,  3.2°S 
dip  latitude)  and  Arecibo,  Puerto  Rico  (18.6'N,  bb.S'W,  30.0°N  dip  latitude).  The  dip 
latitude  is  the  north-south  component  of  the  geomagnetic  coordinate  system.  For 
example,  the  magnetic  dip  latitude  is  90°N  and  90°S  at  the  magnetic  north  and  south 
poles,  respectively.  The  dip  latitude  at  the  magnetic  equator  is  0  degrees.  The  low 
latitude  ionosphere  incorporates  the  unique  neutral  and  electrodynamic  properties  of  the 
equatorial  region  as  well  as  the  dynamics  found  in  the  transition  region  between  the  low 


and  middle  latitudes. 
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The  fundamental  equation  governing  the  amount  of  ionization  at  any  point  in  the 
ionosphere  is  the  continuity  equation: 

^+y.{Nv)  =  P-L  (1) 

at 

Here  N  is  the  electron  density,  V  -{nv)  is  the  transport  term,  P  is  the  production  rate, 
and  L  is  the  loss  rate. 

In  the  low  latitude  ionosphere,  the  E-region  ionization  is  produced  mainly  from 
O2  molecules  by  photoionization: 

O2  +  hv  — >  O2  +  6  (2) 

The  E-region  ion  production  is  balanced  by  the  loss  of  ionization,  principally  by  the 
dissociative  recombination  reactions: 

O;  +  e  ^  O  -h  O  (3) 

and 

NO*  +  e  N  4-  O  (4) 

The  recombination  time  constant  for  the  reactions  (3)  and  (4)  is  3x10"’  cm^sec'* 
[Tascione,  1988].  The  rate  of  change  of  the  E-region  electron  density  is  simply  the 
difference  between  the  ion  production  and  the  ion  loss  processes. 

In  the  F-region  the  principal  ion  is  O*  formed  by 

O  +  /iv  ->  O*  +  e  1027  A  (5) 

and  to  a  lesser  extent 

N2  +  hv  — ^  N2  +  e 


X<196A 


(6) 
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The  loss  process  for  the  F-region  combines  atom-ion  interchange  and  dissociative 
recombination  reactions.  First,  the  following  atom-ion  interchange  reactions  occur: 

+  O2  ^  O2"  +  O  (7) 

and 

O"  +  Nj  ->  NO"  +  N  (8) 

Second,  the  ions  are  lost  by  the  same  dissociated  recombination  processes  outlined  in 
reactions  (3)  and  (4). 

The  low  latitudes  exhibit  unique  transport  effects  due  to  electromagnetic 
dynamics  generated  by  the  thermospheric  winds.  To  a  first-order  approximation  the  F- 
region  neutral  winds  are  driven  by  pressure  gradients  forces  that  arise  from  solar  EUV 
heating  at  F-region  altitudes.  These  pressure  gradients  cause  the  wind  to  blow  from  the 
subsolar  point  to  the  night  side.  Over  Arecibo,  these  winds  and  ambipolar  diffusion 
contribute  to  the  transport  of  ionization  by  moving  the  electrons  and  ions  along  the 
magnetic  field  lines.  During  the  day,  the  pressure  gradients  force  the  winds  in  the 
poleward  direction,  causing  the  ions  to  flow  down  the  magnetic  field  lines.  At  night,  the 
winds  blow  towards  the  equator  forcing  the  ions  to  flow  up  the  field  lines.  These 
nighttime  equatorward  winds  lift  the  level  of  ionization  into  a  region  of  lower  loss  rate 
making  it  possible  to  maintain  the  electron  density  in  the  absence  of  solar  production 
{Walker,  1988]. 

At  equatorial  latitudes,  the  vertical  ExB  drift  is  the  primary  transport  process 
that  affects  the  electron  and  ion  density  distribution.  Data  obtained  by  the  incoherent 
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scatter  radar  at  Jicamarca,  Peru  [e.g.,  Fejer,  1991,  1981;  Fejer  et  al.,  1991;  Kelley, 
1989;  Meriwether  et  al,  1986]  have  provided  some  insight  into  the  effects  of  the 
vertical  drifts  on  the  equatorial  neutral  winds  and  the  electron  density  profile.  The 
effects  of  the  vertical  drifts  on  the  winds  over  the  magnetic  equator  will  be  covered  in  a 
later  section. 


2.3.  UPPER  ATMOSPHERIC  WINDS 

2.3.1.  Stratosphere  and  Mesosphere  Winds 

Hgures  3  and  4  show  the  average  zonal  wind  profile  obtained  from  the  HWM93 
[Hedin  et  al,  1993]  for  the  upper  atmosphere.  These  wind  profiles  are  for  altitudes 
from  the  surface  to  120  km.  Generally  the  winds  are  weakest  at  the  tropopause  (~  12 
km)  and  the  mesopause  (~  85  -  90  km),  and  strongest  at  the  stratopause  (~  55  -  60  km). 

The  general  equation  governing  the  atmospheric  wind  motion  is  an  extension  of 
Newton's  second  law  of  motion  [Hargreaves,  1992]: 

^+2(0x17)  =  — Vp  +  g  +  F  (9) 

dt  p 

The  terms  starting  from  the  left  side  are  the  rate  of  change  of  velocity,  the  Coriolis  term, 
the  pressure  gradient,  acceleration  due  to  gravity,  and  an  acceleration  term  due  to  any 
other  force  that  acts  on  the  medium  [Hargreaves,  1992].  To  a  first  approximation,  the 
stratospheric  and  mesospheric  winds  result  from  the  balance  between  the  Coriolis  force 
and  the  pressure  gradient  force.  This  balance  between  the  Coriolis  force  and  the 
pressure  gradient  force  is  called  the  geostrophic  balance.  The  average  result  for  these 
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Figure  3.  General  magnitude  of  middle  latitude  zonal  winds  from  the  surface  to  120  km 
for  December  at  0000  LT.  The  zonal  wind  data  were  provided  by  the  HWM93 
empirical  model  [Hedin  et  al.,  1993]  for  low  solar  flux  (<I>  =  75)  quiet  geomagnetic 
activity  (Ap  =  4). 
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Figure  4.  Latitude  vs.  altitude  cross-sections  of  the  average  zonal  wind  at  0000  LT. 
The  zonal  wind  data  were  provided  by  the  HWM93  empirical  model  for  low  solar  flux 
(^>  =  75)  and  quiet  geomagnetic  activity  (Ap  =  4).  The  positive  contours  represent  the 
eastward  winds  and  the  shaded  area  or  negative  contours  identify  the  westward  winds. 
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regions  is  a  zonal  wind  —  eastward  in  the  summer,  westward  in  the  winter  [Hargreaves, 
1992],  In  the  stratosphere  and  mesosphere,  the  pressure  gradient  is  related  to  the 
temperature  gradient;  hence,  the  winds  are  sometimes  called  the  thermal  winds. 

Semi-diurnal  and  diurnal  tides  influence  upper  mesospheric  winds.  Solar  EUV 
heating  of  the  ozone  layer  in  the  stratosphere  results  in  tidal  motions  in  both  the  upper 
stratosphere  and  mesosphere.  These  tides  then  propagate  up  through  the  atmosphere, 
influencing  both  the  upper  mesospheric  and  thermospheric  winds. 

2.3.2.  Thermospheric  Winds 

Why  are  the  thermospheric  winds  important  in  earth's  atmosphere?  The 
thermospheric  winds  affect  the  height  of  the  F  layer  and  also  provide  a  source  of  energy 
that  generates  the  electric  fields  [Kelley,  1989].  For  this  thesis,  only  the  thermospheric 
winds  above  200  km  ai’e  studied.  Early  measurements  of  winds  above  200  km  revealed 
a  wind  system  generated  by  pressure  gradients  set  up  by  atmospheric  heating.  This 
heating  was  created  by  absorption  of  solar  EUV  by  the  abundant  species  at  ~  250  km, 
atomic  oxygen.  The  thermal  expansion  of  the  atmosphere  forms  the  “diurnal  bulge”  that 
follows  the  subsolar  point  by  ~  1  -  2  hours  [Rishbeth,  1972].  It  is  the  atmospheric 
expansion  that  gives  rise  to  horizontal  pressure  gradients  that  can  drive  the  horizontal 
wind  system;  hence,  the  primary  forcing  terra  for  thermospheric  wind  patterns  is  the 
pressure  gradient  term.  Although  the  pressure  term  is  the  primaiy  force,  equation  (9) 
indicates  other  forcing  terms  that  sometimes  cannot  be  neglected. 

As  stated  earlier,  the  primary  forcing  term  for  upper  atmospheric  winds  is  the 
pressure  gradient  associated  with  the  diurnal  bulge  that  results  from  solar  EUV  heating. 
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The  diurnal  bulge  is  located  at  low  latitudes  and  is  coincident  with  the  subsolar  point. 
Figure  5  provides  a  first  order  approximation  to  the  thermospheric  winds  at  December, 
March,  and  June  at  250  km  for  moderate  solar  flux  and  quiet  geomagnetic  conditions. 
Overall,  the  average  wind  is  poleward  during  the  day,  eastward  after  sunset,  and 
equatorward  after  midnight. 

A  more  detailed  equation  of  motion  is: 

p— =  -Vp  +  pg  +  TiV^^-V-jc„-2paxf7-pv„,.(f7-l^)  (10) 

dt 

where  p  -  neutral  density,  p  -  pressure,  g  -  acceleration  due  to  gravity,  -  momentum 
flux  density  tensor,  O  -  angular  velocity,  Vni  -  neutral/ion  collision  frequency,  tJ  - 
neutral  velocity,  and  -  ion  velocity.  The  terms  starting  from  the  left  side  are  the  total 
force  on  the  neutral  particles,  the  pressure  gradient  force,  acceleration  due  to  gravity, 
the  viscous  force,  advective  time  derivative  of  the  momentum,  the  Coriolis  force,  and 
the  ion  drag  force.  For  thermospheric  winds,  neglecting  the  Coriolis  force,  the  primary 
limitation  on  the  wind  speed  comes  from  the  various  drag  terms  [Kelley,  1989].  Figure 
6  is  a  graph  showing  the  relative  daytime  magnitudes  of  a  few  terms  in  equation  10. 

Besides  the  pressure  gradient  force,  the  ion  drag  also  affects  the  upper 
thermospheric  winds.  The  ion  drag  becomes  important  when  the  ion  collision 
frequencies  are  much  smaller  than  the  ion  gyrofrequencies.  Examples  of  typical  F- 
region  thermospheric  daytime  ion  collision  frequencies  (Vj)  are  20  s"l  and  ion 
gyrofrequencies  (cq)  are  typically  210  rad/s.  When  the  ion  gyrofrequency  is  greater  than 
the  ion  collision  frequency,  the  ions  gyrate  along  the  magnetic  field  lines;  thus,  the  ions 


Figure  5.  The  temperatures  (solid  curves)  and  winds  (arrows)  at  250  km  for  December, 
March,  and  June.  These  data  were  obtained  using  the  HWM93  for  low  flux  (<E>  =  75) 
and  quiet  geomagnetic  activity  (Ap  =  4) 
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Figure  6.  Graphs  of  neutral  acceleration  time  constants  and  F  for  midday  at  sunspot 
minimum.  The  upper  scale  applies  to  the  acceleration  due  to  the  pressure  gradient  F. 
The  lower  scale  refers  to  the  ion  drag  parameter  v„;;the  Coriolis  parameter  for  latitude 
of  45^,  f,  and  the  normalized  kinematic  viscosity  parameter  Ti/prf  [After  Rishbeth, 
1972.  Reproduced  with  permission  of  Pergamon  Press.] 


cannot  flow  across  the  magnetic  field  lines  like  the  neutral  winds.  In  the  F-region,  the 
effect  of  the  ions  being  confined  to  the  magnetic  field  line  is  a  slowing  of  the  neutral 
winds.  This  “slowing  down”  force  is  called  the  ion  drag. 

A  poleward  wind  forces  the  ions  to  gyrate  down  the  magnetic  field  fines  while  an 
equatorward  wind  causes  the  ions  to  gyrate  up  the  magnetic  field  fines.  If  the  magnetic 
field  fine  is  horizontal  to  the  earth's  surface  (dip  angle  ~  0),  the  winds  will  not  have  any 
effect  on  raising  or  lowering  the  plasma.  This  is  the  situation  at  equatorial  latitudes. 

If  the  plasma  is  lifted  up  the  field  fines,  the  ion  density  is  reduced,  allowing  for 
the  possibility  of  strong  zonal  winds.  Similarly,  if  the  plasma  is  forced  down  the  field 
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lines,  the  ion  density  increases  and  the  neutral  wind  speed  decreases.  Figure  7 
shows  the  effects  of  raising  (lowering)  the  plasma  due  to  equatorward  (poleward) 
winds. 

At  equatorial  latitudes  the  meridional  winds  do  not  raise  or  lower  the  plasma; 
hence,  the  vertical  ion  motion  is  not  affected  by  meridional  winds  at  Arequipa.  The  next 
section  addresses  the  wind  dynamics  at  the  equatorial  latitudes. 

2.3.3.  Equatorial  Thermospheric  Winds 

The  last  section  examined  the  effect  of  the  steady  drag  on  the  neutral  wind  for 
the  case  where  the  ion  gyrofrequency  is  greater  than  the  ion  collision  frequency.  At 
equatorial  latitudes  the  meridional  effect  on  raising  and  lowering  the  ions  along  the  field 
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Figure  7.  Illustration  of  how  the  thermospheric  winds  affect  the  ionosphere. 
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lines  is  negligible  since  the  magnetic  field  is  nearly  parallel  to  the  earth's  surface  (dip 
angle  =  0).  Before  exploring  further  the  wind  dynamics  of  the  equatorial  latitude 
region,  a  short  discussion  is  presented  on  the  theoretical  wind  speed  calculations  at  this 
region. 

At  equatorial  latitudes  the  Coriolis  force  is  neglected  but  the  drag  force  is  still 
important.  Thus,  for  steady-state  conditions,  the  equation  of  motion  (10)  reduces  to 
[Kelley,  1989]: 

^=F-v,ft=Q  (11) 

dt 

where  F  is  the  pressure  gradient  force  and  v„,t7  is  the  drag  force.  Equation  (11) 
indicates  that  the  equation  of  motion  reduces  to  a  balance  between  the  pressure  gradient 
and  the  drag  force.  Solving  for  the  velocity,  the  neutral  velocity  becomes: 

17=-—  (12) 

V„i 

By  using  the  values  in  Figure  6,  the  pressure  gradient  and  viscosity  coefficient  can  be 
obtained  as  a  function  of  altitude.  For  example,  at  300  km,  the  velocity  of  the  neutral 
wind  is  calculated  to  approximately  130  m/s.  This  result  is  well  below  early  calculations 
of  the  neutral  wind  by  Lindzen  [Kelley,  1989].  Because  Lindzen  did  not  include  the 
drag  force,  his  calculated  wind  speed  was  550  m/s  compared  to  the  more  realistic  130 
m/s.  Thus  the  ion  drag  is  the  crucial  factor  that  provides  correct  results  for  typical 
upper  thermospheric  wind  speeds. 

What  causes  the  strong  eastward  zonal  winds  near  the  magnetic  equator  during 
the  hours  of  2000-2400  LT?  Without  the  effect  of  meridional  winds  on  raising  or 
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lowering  the  ions,  some  other  mechanism  is  necessary  to  explain  the  upper 
thermospheric  wind  pattern  during  the  postsunset  hours.  The  prereversal  enhancement 
of  the  zonal  electiic  field  is  indirectly  responsible  for  the  strong  easterly  winds  [Kelley, 
1989].  The  effect  of  the  prereversal  enhancement  is  to  introduce  a  strong  vertical 
ExB  that  lifts  the  ions  during  1700-2000  UT  (see  Figures  8  and  9).  Figure  9  shows 
the  upward  vertical  drift  is  most  pronounced  during  the  equinox  and  June  solstice 
periods  and  has  large  day-to-day  and  solar  cycle  vertical  drift  shortly  before  midnight, 
which  increases  the  ion  drag  and  reduces  the  variations  [Fejer  et  al,  1991].  With  the 
ions  out  of  the  way,  the  ion  drag  is  greatly  reduced,  thus  nearly  doubling  the 
thermospheric  wind  speeds  (see  Figure  10).  This  effect  is  clearly  shown  in  the  zonal 
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Figure  8.  Calculated  vertical  ion  drift  velocities  for  several  driving  wind  components. 
The  soUd  line  includes  only  the  tidal-driven  E-region  dynamo,  while  the  dashed  line 
includes  the  F-region  dynamo  as  weU.  Typical  vertical  plasma  drifts  are  indicated  by  the 
dotted  line  [  After  Heelis  et  al,  1974.  Reproduced  with  permission  of  Pergamon  Press.] 
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Figure  9.  Measured  vertical  (positive  for  upward  motions)  plasma  drifts  in  the 
equatorial  F-region  for  various  seasons  and  solar  cycle  conditions. 
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Figure  10.  Contours  of  the  zonal  neutral  wind  at  250  km,  calculated  from  the 
thermospheric  model  by  Anderson  and  Roble  (1974)  as  a  function  of  latitude  and  local 
time;  (a)  without  the  postsunset  enhancement  of  the  ExB  drift  velocity  and  (b) 
including  the  postsunset  enhancement  [After  Anderson  and  Roble,  1974.  Reproduced 
with  permission  of  the  American  Geophysical  Union]. 
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wind  measurements  made  by  the  Dynamics  Explorer  2  satellite  (see  Figure  11). 

The  previous  two  sections  provided  only  a  glimpse  of  how  the  upper 
thermospheric  winds  are  generated  and  effected  by  ion  drag  and  ExB  vertical  drifts. 
Other  ionospheric  effects  such  as  gravity  waves,  tides,  and  geomagnetic  storms  also 
play  a  significant  role  in  the  dynamics  of  the  upper  thermosphere.  While  this  work  does 
not  highlight  the  effects  of  gravity  waves  and  tides,  it  is  important  to  realize  that  the 
vertical  drifts  and  ion  drag  are  not  the  only  mechanisms  affecting  the  thermospheric 
neutral  winds.  Chapter  6  will  discuss  the  effects  of  geomagnetic  storms,  also  known  as 
dismrbed  conditions,  on  the  low  latitude  neutral  winds. 

2.4.  PREVIOUS  WIND  STUDIES 

Several  climotological  studies  of  the  thermospheric  winds  have  been  conducted 
at  Arequipa  and  Arecibo.  Although  significant  longitudinal  variations  of  the  neutral 
winds  were  noted  in  the  Dynamic  Explorer  2  satellite  measurements  [Raghavarao  et  al, 
1991],  the  measured  neutral  winds  over  Arequipa  and  Arecibo  were  shown  to  be 
consistent  with  other  low  latitude  sites  at  different  longitudes  [Biondi  et  al,  1991,  1990; 
Meriwether  et  al,  1986;  Sahai  et  al,  1992].  This  section  briefly  presents  the  most 
significant  results  of  these  studies.  We  will  also  highlight  the  magnetic  activity  smdies 
since  this  work  focus  on  the  neutral  wind  response  to  storm  time  events. 

For  Arequipa,  results  by  Biondi  et  al  [1990,  1991]  and  Meriwether  et  al 
[1986]  showed  the  thermospheric  winds  exhibited  both  solar  cycle  effects  and 
seasonal  variations,  with  the  seasonal  variations  being  more  pronounced  (see  Figures  12 
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Figure  11.  Dynamics  Explorer  2  data  of  electron  density  (upper  panel),  zonal  winds 
(middle  panel),  and  temperature  (lower  panel)  near  300  km  at  low  latitudes.  The 
measurements  were  taken  on  November  20,  1981  at  19(X)  hrs  LST  (left  graph)  and 
November  20,  1982  at  18.9  hrs  LST  (right  graph).  These  graphs  show  an  inverse 
relationship  between  the  magnitudes  of  the  wind  velocities  (‘Z’)  and  the  electron  density 
(‘P’)  profiles.  Note  at  the  Dip  Equator,  the  measurements  support  the  increased  zonal 
winds  at  the  time  of  the  sharp  increase  of  the  upward  velocity  in  the  dusk  sector.  [After 
Raghavarao  et  al,  1991.  Reproduced  with  permission  of  the  American  Geophysical 
Union.] 
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and  13).  These  studies  did  not  address  magnetic  activity  effects. 

Thermospheric  wind  studies  over  Arecibo  by  Burnside  and  Tepley  [1989]  and 
Burnside  et  al.  [1981,  1983,  1987,  1991]  showed  the  neutral  winds  are  not  affected  by 
changes  in  the  solar  cycle  activity,  although  there  exist  pronounced  seasonal  variations. 
This  study  did  not  provide  a  detailed  analysis  of  the  solar  cycle  effects  on  the  seasonal 
averages  (see  Figure  14).  As  for  storm  time  effects,  these  studies  concentrated  on  a  few 
disturbed  days  and  did  not  address  the  issue  of  the  repeatability  of  the  neutral  wind 
response.  During  a  geomagnetic  disturbance  on  the  night  of  14  January  1988  (Kp  >  7), 
Burnside  et  al.  [1991]  showed  measured  FPI  winds  over  Arecibo  (see  Figure  15).  The 
results  in  Figure  15  show  mainly  westward  neutral  winds  and  near  zero  meridional 
winds  during  the  night.  These  results  also  show  a  return  to  average  eastward  zonal  and 
poleward  meridional  wind  during  the  night  of  15  January  where  the  three  hourly  Kp 
returned  to  moderate  to  quiet  geomagnetic  conditions  (Kp  <  3). 
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Figure  12.  Arequipa  equinox  month  average  meridional  and  zonal  winds  from  1983  to 
1990.  The  short-dash  curves  labeled  H  represent  predictions  of  the  HWM-87  by  Hedin 
et  al.  [1988].  [After  Biondi  et  al.,  1991,  Reproduced  with  permission  from  the 
American  Geophysical  Union]. 
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Figure  13.  Arequipa  June  solstice  average  meridional  and  zonal  winds  from  1983  to 
1990.  The  months  included  in  each  year’s  average  are  June  and  July,  except  for  1984, 
where  May  and  August  are  used,  as  indicated  by  the  M  and  A  symbols.  The  short-dash 
curves  labeled  H  represent  predictions  of  the  HWM-87  by  Hedin  et  al  [1988].  The 
long-dash  curves  labeled  V  are  the  VSH  TGCM  predictions  by  Killeen  et  al.  [1987]. 
[After  Biondi  et  al,  1991,  Reproduced  with  permission  from  the  American  Geophysical 
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Figure  14.  The  annual  variation  of  the  average  northward  meridional  and  eastward 
zonal  Fabry-Perot  winds  over  Arecibo.  The  open  circles  represent  the  results  obtained 
by  sorting  the  measurements  from  1980  to  1983  into  30-day  periods.  The  crosses  show 
the  results  for  the  measurements  from  1984  to  1988.  [After  Burnside  and  Tepley, 
1989.  Reproduced  with  permission  by  the  American  Geophysical  Union]. 
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Figure  15.  Zonal  neutral  wind  velocities  measured  by  the  Fabry-Perot  interferometer 
(circles)  on  four  consecutive  nights  in  January  1988.  Cross  show  the  ion  velocity  at  300 
km  measured  by  the  incoherent  scatter  radar.  The  dashed  lines  are  the  zonal  winds 
given  by  the  HWM-87  empirical  model.  [After  Burnside  et  aL,  1991.  Reproduced  with 
permission  by  the  American  Geophysical  Union]. 
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CHAPTER  3 

MODELING  TECHNIQUES 

3.1.  INTRODUCTION 

The  long-term  average  and  the  day-to-day  variability  of  the  ionosphere  can  be 
analyzed  using  upper  atmospheric  models.  These  models  further  our  understanding  of 
the  physical  mechanisms  that  contribute  to  the  long-term  and  daily  variability  of  the 
ionosphere.  This  thesis  has  already  presented  several  model  predictions  from  the 
empirical  models  HWM93  [Hedin  et  al.,  1993]  and  MSISE90  [Hedin,  1991],  and  the 
first  principal  model  TIGCM93  [Fesen  et  al,  1993].  This  chapter  provides  the 
historical  aspects  of  upper  atmospheric  modeling  and  outlines  several  models  used  for 
this  study.  The  chapter  concludes  with  a  more  detailed  description  of  the  TIGCM93 
physical  equations  and  low  latitude  neutral  wind  predictions. 

3.2.  UPPER  ATMOSPHERIC  MODELING 

3.2.1.  History  of  Atmospheric  Modeling 

The  development  of  upper  atmospheric  models  (above  20  km  altitude)  is  based 
on  tropospheric  forecast  models  such  as  the  Limited  Fine  Mesh  (LFM)  and  the  Nested 
Grid  Model  (NGM).  The  National  Weather  Service  uses  the  LFM  and  NGM  for  daily 
weather  prediction.  Because  of  the  similarities  in  atmospheric  dynamics  in  both  the 
upper  and  lower  atmosphere,  upper  atmospheric  scientists  borrowed  many  ideas  from 
the  lower  atmospheric  forecast  models  and  incorporated  them  into  upper  atmospheric 
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models.  This  section  provides  a  brief  history  of  lower  atmospheric  modeling  because  of 
its  importance  to  the  thermospheric  modeling  efforts. 

The  British  scientist  L.F.  Richardson  was  the  first  to  attempt  to  predict  weather 
numerically  [Holton,  1992].  His  work  is  considered  the  cornerstone  of  weather 
prediction  and  is  based  on  finite  difference  calculations.  Unfortunately,  Richardson  was 
hampered  by  two  major  problems.  Since  his  work  was  done  in  the  1920s,  he  did  not 
have  the  luxury  of  computers  to  perform  the  numerous  arithmetic  operations. 
Richardson  estimated  that  64,000  people  would  be  required  to  keep  up  with  these 
global  forecasts  [Holton,  1992]. 

Richardson's  second  problem  concerned  his  initial  attempt  to  predict  weather 
between  two  grid  points.  This  first  weather  prediction  failed  so  miserably  that  people 
doubted  that  weather  prediction  would  become  a  feasible  operation.  It  was  not  until 
after  World  War  II  that  researchers  realized  the  problem  associated  with  Richardson's 
first  forecast  attempt.  When  Richardson  set  up  the  differential  equations  for  his 
forecast,  he  included  terms  for  high  speed  sound  and  gravity  wave  solutions.  As  the 
numerical  prediction  was  marched  forward  in  time,  the  sound  and  gravity  wave 
solutions  amplified  so  much  that  they  obliterated  the  correct  numerical  predictions! 

After  1950,  meteorologists  made  significant  progress  in  weather  prediction  and 
modeling.  It  was  not  until  20  years  later,  however,  that  upper  atmospheric  scientists 
began  exploring  the  possibility  of  applying  these  lower  atmospheric  equations  to  the 
regions  above  the  troposphere. 
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In  1979,  the  Theimospheric  General  Circulation  Model  (TGCM)  was  created  to 
better  understand  the  processes  affecting  the  upper  atmospheric  weather.  The  TGCM, 
based  on  similar  equations  used  to  predict  tropospheric  weather,  has  undergone  several 
model  revisions.  These  changes  will  be  highhghted  later  in  this  paper. 

3.2.2.  Upper  Atmospheric  Models 

Several  models  used  to  study  the  upper  atmosphere  include  the  Mass 
Spectrometer  and  Incoherent  Scatter  Extended  (MSISE90)  Model,  the  Horizontal  Wind 
Model  (HWM93),  the  Vector  Spherical  Harmonic  (VSH)  model,  and  the  TIGCM. 
Before  describing  these  models,  the  distinction  should  be  made  between  empirical,  semi- 
empirical,  and  first  principle  models.  An  empirical  model  is  a  model  whose  calculations 
are  based  on  collected  data  and  are  statistical  by  nature.  Examples  of  an  empirical 
model  are  the  MSISE90  and  the  HWM93  since  their  calculations  are  based  on  data 
collected  from  satellite  mass  spectrometer,  incoherent  scatter  radar  measurements  and 
satellite  measurements.  A  physical  model  is  a  model  whose  calculations  are  based  on 
physical  principles.  The  TIGCM  is  an  example  of  a  physical  model.  A  semi-empirical 
model  is  a  model  whose  calculations  are  based  on  collected  data  and  physical  principles. 
The  VSH  model  is  an  example  of  a  semi-empirical  model. 

3.2.2. 1.  Mass  Spectrometer  and  Incoherent  Scatter 
(MSISE90)  Model 

The  MSISE90  is  an  empirical  model  based  on  in  situ  and  Incoherent  Scatter 
Radar  measurements.  The  MSISE90  version  provides  modeled  data  from  the  surface 
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to  approximately  3000  km  altitude.  The  model  calculates  profiles  of  number  densities, 
total  mass  density,  and  neutral  temperature  for  N,  N^,  O^,  O,  He,  Ar,  and  H  [Hedin, 
1991].  These  predictions  are  based  on  data  collected  from  several  satellites,  rocket 
probes  (in  situ),  and  five  incoherent  scatter  radars  [Hedin,  1991].  The  MSISE90  model 
predictions  are  functions  of  several  geophysical  values  such  as  10.7 -cm  solar  flux,  Ap 
index,  latitude  and  longitude,  and  universal  time. 

3.2.2.2.  Horizontal  Wind  Model  (HWM93) 

After  discussing  the  MSISE90  model,  it  is  only  natural  to  wonder  if  a  similar 
type  of  empirical  model  is  used  to  estimate  upper  atmospheric  winds.  The  HWM93  is 
the  MSISE90  equivalent  for  upper  atmospheric  winds.  The  HWM93  is  based  on 
thermospheric  wind  data  obtained  from  the  Atmosphere  Explorer  E  (AE-E)  and 
Dynamics  Explorer-2  (DE-2)  satellites,  and  ground-based  data  from  several  incoherent 
scatter  radars  and  Fabry-Perot  optical  interferometers  [Hedin  et  al,  1991]. 

The  HWM93  model  requires  the  user  to  input  the  day  of  year,  time  of  day  (UT), 
altitude,  geographic  latitude  and  longitude,  local  solar  time,  3-month  average  and 
previous  day  value  of  the  10.7-cm  solar  flux  index,  and  either  the  daily  Ap  magnetic 
index  or  prescribed  history  of  3-hour  Ap  indices  [Hedin  et  al.,  1991].  Once  the  user 
enters  these  data,  HWM93  gives  the  estimated  zonal  and  meridional  winds  from  100  km 
to  500  km. 

T  -ike  the  MSISE90  model,  the  HWM93  uses  two  different  fitting  routines  above 
and  below  200  km.  Above  200  km,  where  viscosity  limits  neutral  wind  latitudinal 
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gradients,  these  altitude  variations  are  each  represented  by  an  extension  of  the  Bates 
formula.  This  is  the  same  formula  used  for  thermospheric  temperature  profiles  [Hedin 
et  al,  1991].  Below  200  km  the  wind  profiles  are  represented  by  a  cubic  spline.  This 
cubic  spline  is  defined  by  cubic  polynomials  between  specified  nodes  with  first  and 
second  derivatives  continuous  across  the  interior  nodes  [Hedin  et  al.,  1991].  These 
nodes  are  located  between  250  km  and  300  km  because  of  the  limited  data  in  this 
altitude  range. 

3.2.2.3.  Vector  Spherical  Harmonic  (VSH)  Wind  Model 

The  VSH  model  is  a  computer  subroutine  that  provides  a  description  of  the 
composition  and  dynamics  of  the  thermosphere  [Killeen  et  al.,  1987].  Given  the  time, 
location,  the  Ap  index,  and  the  10.7-cm  solar  flux,  the  VSH  routine  returns  a  variety  of 
parameters.  First,  the  VSH  routine  returns  densities  for  O,  O^,  N^,  H,  O  ,  electron,  and 
total  mass.  Second,  the  routine  returns  the  mixing  ratios  for  O,  Nj,  and  O2.  The 
dynamic  parameters  returned  are  the  neutral  zonal  winds,  meridional  winds,  vertical 
winds,  neutral  and  ion  temperatures,  and  pressure.  These  data  are  provided  for  altitudes 
ranging  from  130  km  to  600  km  [Killeen  etal,  1987]. 

Two  types  of  data  are  used  to  drive  the  VSH  model:  coefficient  library  and  real¬ 
time  data  [Killeen  et  al,  1987].  The  coefficient  library  consists  of  TTGCM  data  output. 
The  real-time  data  are  made  up  of  solar  flux  and  magnetic  activity  indices,  as  well  as  any 
compositions,  density,  and/or  wind  information  as  obtained  by  satellite  measurement 
[Killeen  et  al,  1987].  The  VSH  model  is  a  semi-empirical  model  because  of  the 
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dependence  on  collected  data  and  the  TIGCM  calculated  data.  Because  the  upper  limit 
of  the  TIGCM  is  500  km,  the  VSH  model  uses  the  MSIS-86  model  results  above  500 
km  altitude. 

Two  reasons  make  the  VSH  model  an  appealing  option.  First,  the  model  makes 
temperature  results  available  for  comparison  with  other  empirical  models  [Killeen  et  al, 
1987].  Second,  the  VSH  model  uses  more  coefficients  in  its  analysis.  This  provides 
more  accurate  results  than  other  empirical  models  such  as  MSIS  [Killeen  et  ah,  1987]. 
The  main  drawback  is  that  this  model  does  not  improve  on  the  TIGCM  results. 

3. 2.2.4.  Thermosphere/Ionosphere  General  Circulation 
Model  (TIGCM93) 

This  following  section  describes  NCAR's  TIGCM93.  To  develop  an 
appreciation  for  this  model,  the  next  section  will  focus  on  the  history,  description, 
calculated  results,  primitive  equations,  and  model  output  of  the  TIGCM93.  For  the 
remainder  of  this  chapter,  we  shall  refer  to  this  model  as  the  TIGCM  rather  than  the 
TTGCM93.  The  “93”  refers  to  the  version  number  of  the  most  current  TTGCM 

3.3.  THERMOSPHERE/IONOSPHERE  GENERAL  CIRCULATION 

MODEL  CnGCM) 


3.3.1.  TTGCM  History 

The  first  NCAR  General  Circulation  Model  was  operational  in  1979.  Dr. 
Raymond  Roble  is  considered  the  father  of  the  thermospheric  general  circulation  models 
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since  he  was  the  first  to  adjust  Holton's  primitive  equations  and  apply  the  revised 
equations  to  thermospheric  physics.  TABLE  2  lists  a  chronological  development  of  the 
thermospheric  general  circulation  models. 

The  first  thermospheric  general  circulation  model,  released  in  1979,  was  called 
the  Thermosphere  General  Circulation  Model  (TGCM).  This  model  used  the  revised 
primitive  equations,  described  in  a  following  section,  to  predict  the  general  circulation 
of  the  thermosphere.  Besides  the  primitive  equations,  the  TGCM  relied  on  several 
empirical  models  to  supplement  its  calculations.  The  empirical  models  (see  TABLE  2) 
provided  input  for  the  dynamo  processes  for  both  E-  and  F-regions,  neutral  gas 


TABLE  2.  Chronological  Development  of  the  TGCM. 


Year 

Empirical  Model 

General  Circulation  Model 

1979 

Ionospheric 

TGCM 

Dynamo 

MSIS 

Aurora 

Solar 

Tides 

(100 -500  km) 

1989 

Dynamo 

TIGCM 

Aurora 

Solar 

Tides 

(95 -500  km) 

1991 

Aurora 

TIE-GCM 

Solar 

Tides 

(95  -  500  km) 

1992 

Aurora 

TIME-GCM 

Solar 

(30 -500  km) 
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concentration  (MSIS),  ionospheric  coupling  interactions,  and  external  forcing 
mechanisms  (aurora,  solar,  and  tidal  motions). 

The  TIGCM  was  released  in  1989.  The  differences  between  the  TIGCM  and  the 
TGCM  were  the  inclusion  of  coupling  mechanisms  in  the  ionosphere  and  prediction  of 
neutral  concentrations.  Since  the  TIGCM  predicted  the  ionospheric  coupling  and 
neutral  gas  concentrations,  the  requirement  for  the  Ionospheric  and  MSIS  empirical 
models  was  dropped.  In  addition,  the  lower  boundary  of  the  model  was  dropped  to  95 
km.  However,  this  model  has  not  been  fully  tested. 

3.3.2.  TIGCM  Description 

The  TIGCM  model  treats  all  three  spatial  dimensions  (x,  y,  and  z),  including 
time,  by  finite  differences,  allowing  more  flexibility  than  thermospheric  models  prior  to 
the  TGCM  [Dickinson  et  ai,  1981].  The  TIGCM  solves  the  primitive  equations  of 
dynamic  meteorology  but  adapts  to  the  physics  appropriate  to  thermospheric  heights. 

The  TIGCM  contains  most  of  the  physics  to  model  a  coupled  ionosphere  and 
thermosphere  atmosphere.  Besides  the  primitive  equations  (thermodynamic,  east  and 
west  momentum,  continuity,  and  hydrostatic  equations),  the  model  includes  most  of  the 
basic  chemical  reactions  that  occur  in  the  thermosphere.  The  most  apparent  limitation 
of  this  model  is  that  it  does  not  contain  the  F-region  dynamo  and  as  a  result  cannot  be 
considered  a  true  physical  model. 

The  TIGCM  calculates  global  distributions  of  the  neutral  gas  temperatures; 
zonal,  meridional,  and  vertical  winds;  the  major  neutral  gas  constituents  number 
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densities  of  O2,  O,  and  N^;  and  the  minor  neutral  gas  constituent  number  densities  of 
N(^D),  NC'^S),  and  NO  [Roble,  1992].  The  interactive  model  of  the  ionosphere  solves 
for  global  distributions  of  the  ion  species  O^,  Oj,  Nj,  NO^,  and  the  electron 
density.  Additional  predictions  for  the  interactive  ionosphere  model  include  electron 
and  ion  temperatures.  The  model  has  an  effective  5x5  latitude  by  longitude  horizontal 
grid  with  25  constant  pressure  surfaces  extending  approximately  between  95  and  500 
km  in  altitude  with  a  vertical  resolution  of  two  grid  points  per  scale  height  [Roble, 
1992].  The  model  time  step  is  5  minutes  and  it  takes  about  20  minutes  for  the  NCAR 
CRAY  CPU  to  simulate  1  day  [Roble,  1992]. 

The  TIGCM,  tike  the  TGCM,  uses  several  empirical  models  for  its  calculations. 
These  empirical  models  are  listed  in  TABLE  3.  If  no  name  is  specified  for  the  empirical 
model,  the  author’s  name  is  used  (italics);  otherwise  the  name  of  the  model  is  listed. 
TABLE  3  indicates  that  the  TIGCM  model  is  dependent  on  empirical  models  to 


TABLE  3.  Empirical  Models  Used  in  the  TIGCM 


Empirical  Model 

Model  Calculations 

Richmond 

low  and  middle  latitude  electric  fields  [Richmond  et  al., 
1980] 

Heelis 

high  latitude  electric  fields  [Heelis  et  al.,  1982] 

FuUer-Rowell  and  Evans 

auroral  particle  inputs  [Fuller-Rowell  and  Evans,  1987] 

Hinteregger 

solar  EUV  inputs  [Hinteregger,  1981] 

Forbes  and  Vial 

amplitudes  and  phases  of  the  upward  propagating 
semidiurnal  tides  [Forbes  and  Vial,  1989, 1991] 
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successfully  predict  upper  atmospheric  properties. 


3.3.3.  TIGCM  Model  Equations 

The  equations  used  by  the  TIGCM  are  based  on  the  pressure  coordinate 
primitive  equations  of  lower-atmospheric  dynamic  meteorology  [Dickinson  et  ai, 
1981].  The  primitive  equations  are  the  thermodynamic,  eastward  and  northward 
momentum,  continuity,  and  hydrostatic  equations. 


3.3.3.I.  Thermodynamic  Equation 

The  thermodynamic  equation  in  spherical  coordinates  is 
dT_  ge^  d  (K^dr 


dt  p„c.  dz 


H  dz. 


-aT 


-VVT-w 


dT  RT 


S  + 


V 


dz  CpHlj 


Q' 


(13) 


where  T  -  temperature,  z  -  loggCpo/p).  S  "  acceleration  of  gravity,  -  reference 
pressure  50pPa,  Cp  -  specific  heat  per  unit  mass,  Kj.  -  thermal  conductivity  coefficient, 
H  -  mean  scale  height,  a  -  radiative  damping  coefficient,  V  •  V  -  horizontal  advective 
operator,  w  -  vertical  motion,  S  -  static  stability,  m  -  specified  mean  molar  mass,  and  Q' 
-  heating  term.  The  terms  from  left  to  right  in  the  equation  are  the  total  heating  rate, 
downward  heat  conduction,  radiative  cooling,  heat  advection,  the  adiabatic  heating  and 
cooling  due  to  expansion  and  contraction  of  the  gases  in  the  thermosphere,  and  the 
heating  term.  The  heating  term  includes  solar  EUV  and  UV  heating.  Joule  heating,  and 


41 

heating  due  to  auroral  particle  precipitation.  Mathematically,  this  Q'  term  is: 

— =-(a+fi;+e,)  (w) 

where  -  rate  of  solar  heating  per  unit  mass,  Qj  -  rate  of  Joule  heating  per  vmit  mass, 
and  Qp  -  rate  of  auroral  particle  precipitation  heating  per  unit  mass 


3.3.3.2.  Momentum  Equations 

The  TIGCM  uses  the  eastward  and  northward  momentum  equations  in  spherical 


coordinates  are: 

Eastward  Momentum  Equation 


du  _ 

ge^  d  f 

Tt~ 

p.  du 


+|-yVM-— tan(t) 


+  fv  +  ['>^„  (Ml  -u)  +  X^  (vj  -  v)} 

1  1  ^ 

I  rcos(t)  dX  dz 


(15) 


Northward  Momentum  Equation 

A 


dv  ge^  d  f  [L  du 


dt  Pq  dz 


\H  dzj 
2 


-fU+{Xyy{Vj  -V)+X^^(H-M^)} 


(16) 


u  I  1 B*!*'  dv 

-V  •  Vv-  —  tan  4)  - -w — 

r  1  r  3<1)  dz 


where  u  -  eastward  neutral  velocity,  z  -  loge(po/p),  v  -  northward  neutral  velocity,  p,  - 
coefficient  of  viscosity,  g  -  acceleration  of  gravity,  -  reference  pressure  50p.Pa,  Uj  - 
zonal  ion-drift  velocity,  X^,X^  ^'Xy^,Xy^  -  ion  drag  tensor,  H  -  mean  scale  height,  V;  - 
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meridional  ion-drift  velocity,  y  •  V  -  horizontal  advective  operator,  w  -  vertical  motion, 
<I>  -  geopotential,  f  -  Coriolis  parameter,  and  (j)  -  latitude. 

The  terms  from  left  to  right  (for  both  momentum  equations)  are  the  total  net 
force  per  unit  mass,  force  per  unit  mass  due  to  vertical  viscosity,  Coriolis  for  per  unit 
mass,  ion-drag  force  per  unit  mass,  nonlinear  momentum  advection  force  per  unit  mass, 
and  the  pressure  force  per  unit  mass. 

3.3.3.3.  Continuity  Equation 

The  continuity  equation  in  spherical  coordinates  is 

— ^  ,  .f.  (vcos(l))+ — =  0 
rcos(t)a<t)  rcos<poA.  az 

where  <[)  -  latitude,  r-radial  distance  from  the  earth,  v  -  northward  velocity,  u  -  eastward 

velocity,  w  -  dz/dt  (Vertical  motion'),  z  -  log^—  vertical  height  in  pressure  coordinates, 

P 

and  "k  -  longitude. 

The  continuity  equation  expresses  the  conservation  of  mass  for  a  fluid. 

3.3.3.4.  Hydrostatic  Equation 
The  hydrostatic  equation  is: 

^  =  R{T^+T)lm  (18) 

az 

where  <!>'  -  geopotential,  R  -  universal  gas  constant,  Tq  -  global  mean  temperature,  T  - 


temperature,  and  m  -  mass. 
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The  ideal  gas  law  and  the  hydrostatic  balance  combine  to  form  the  hydrostatic 
equation.  The  ideal  gas  law  for  the  atmosphere  describes  the  thermodynamic  state  of 
the  atmosphere  at  any  point  and  is  determined  by  the  values  of  pressure,  temperature, 
and  density  at  a  specified  point  [Holton,  1992]: 

p  =  pRT  (19) 

where  p  -  pressure,  p  -  atmospheric  density,  R  -  gas  constant,  and  T  -  temperature. 

The  hydrostatic  balance  condition  is  a  good  approximation  for  the  vertical 
dependence  of  the  pressure  field  in  the  real  atmosphere.  The  hydrostatic  balance 
condition  defines  the  relationship  between  pressure  and  height 

There  is  a  difference  between  geopotential  and  geopotential  height.  The  gradient 
of  a  potential  function  O,  the  geopotential,  is  equal  to  the  acceleration  due  to  gravity 
[Holton,  1992]: 

VO  =  -g  (20) 

By  setting  the  geopotential  to  zero  at  mean  sea  level,  the  geopotential  at  height  z 
becomes  the  work  required  to  raise  a  unit  mass  to  height  z  from  mean  sea  level  {Holton, 
1992): 


0 

For  most  atmospheric  calculations,  the  geopotential  height 


(21) 


is  used  instead  of  the 


geopotential  <E>.  The  geopotential  height  is  mathematically  defined  as 


go 


3.3.4.  TIGCM  Diagnostic  Processor 

One  of  the  outputs  of  the  NCAR  TIGCM  is  a  history  file  containing  the 
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calculated  wind,  temperature,  and  composition  fields  at  all  model  grid  points  for  each 
hour  [Killeen  and  Roble,  1984].  The  TIGCM  diagnostic  processor  is  run  after  the 
TIGCM  completes  its  calculations.  The  processor  reads  the  history  file,  ion  drag 
tensors,  ion  drift  velocities,  specific  heats,  coefficients  of  viscosity,  and  coefficients  of 
thermal  conductivity,  and  then  applies  horizontal  and  vertical  finite  differencing  to 
obtain  the  individual  forcing  terms  appearing  in  the  momentum  equations  at  a  given 
altitude  for  a  given  time  [Killeen  and  Roble,  1984].  The  processor  provides  these  terms 
both  for  the  global  TIGCM  grid  at  the  given  universal  time  and  for  selected  grid  points 
for  each  hour  of  local  time  [Killeen  and  Roble,  1984]. 

3.4.  TIGCM  PREDICTIONS 

Figures  16  and  17  show  the  TIGCM  predictions  for  the  neutral  winds  over 
Jicamarca,  Peru  (17.5”S,  71.0“W,  dip  latitude  ~  1.3'’S).  Because  the  magnetic  field  lines 
are  horizontal  over  Jicamarca,  the  effects  of  the  meridional  winds  on  the  motion  of  the 
ions  are  negligible.  Section  2.3.2  explained  how  the  prereversal  enhancement  or  the 
vertical  ExB  drift  drives  the  plasma  upward  thereby  reducing  the  ion  drag.  We  note 
that  the  TIGCM  predicted  winds  are  consistent  with  this  theory  since  we  find  the 
maximum  winds  occur  during  2000  LT  -  2300  LT. 

Figures  18  and  19  show  the  output  for  the  TIGCM93  meridional  and  zonal 
winds  for  Arecibo,  Puerto  Rico  (18°N,  71”W).  The  figures  show  the  December  and 
June  solstice  winds  at  300  km  for  solar  minimum  conditions  (O  =  75).  The  most 
noticeable  effect  is  the  north-to-south  wind  shift  between  the  December  and  June 
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JICAMARCA  Northward  Winds  (TIGCM93) 
O  200km  A  300km  +  400km 


Figure  16.  Meridional  winds  from  TIGCM93  model  for  Solar  Minimum  Conditions. 
F10.7  =  75,  Northward  is  positive,  and  local  time  is  centered  about  2400  LT.  Model 
results  for  grid  point  near  Jicamarca  (17.5"S,  70.0°W)  for  December  Solstice. 


JICAMARCA  Eastward  Winds  (TIGCM93) 
O  200km  A  300km  +  400km 


Figure  17.  Zonal  winds  from  TIGCM93  model  for  Solar  Minimum  Conditions.  F10.7  = 
75,  Eastward  is  positive,  and  local  time  is  centered  about  2400  LT.  Model  results  for 
grid  point  near  Jicamarca  (17.5°S,  70.0”E)  for  December  Solstice. 
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ARECIBO  Northward  Winds  at  250km  (TIGCM93) 
O  December  Solstice  A  June  Solstice 


Figure  18.  Meridional  winds  from  TIGCM93  model  for  Solar  Minimum  Conditions. 
F10.7  =  75,  Northward  is  positive,  and  local  time  is  centered  about  2400  LT.  Model 
results  for  grid  point  near  Arecibo  (17.5°N,  bS.O'W)  for  both  December  and  June 
Solstice. 


ARECIBO  Eastward  Winds  at  250km  (TIGCM93) 
O  December  Solstice  A  June  Solstice 


Figure  19.  Zonal  winds  from  TIGCM93  model  for  Solar  Minimum  Conditions.  F10.7 
=  75,  Eastward  is  positive,  and  local  time  is  centered  about  2400  LT.  Model  results  for 
grid  point  near  Arecibo  (17.5°N,  65.0"W)  for  both  December  and  June  Solstice. 
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solstices.  This  wind  change  is  a  result  of  the  sun's  apparent  path  across  the  earth's 
surface.  For  the  December  solstice,  the  sun's  position  crosses  south  of  Arecibo  while 
during  the  June  solstice  the  sun's  position  crosses  north  of  Arecibo.  The  pressure 
gradients  over  Arecibo  are  effected  by  the  sun's  apparent  path;  hence,  these  pressure 
gradients  are  directed  in  such  a  way  as  to  cause  predominate  wind  situations  for  the 
December  and  June  solstice  plots.  The  main  differences  between  the  modeled  zonal 
winds  are  the  smaller  eastward  wind  speeds  during  the  night  for  the  June  solstice 
conditions  and  earlier  morning  reversal  time  during  the  summer. 

Considering  only  the  ion  drag  and  viscosity  forces,  a  poleward  wind  forces  the 
plasma  to  move  down  the  field  lines  while  the  equatorward  wind  forces  the  plasma  to 
move  up  the  field  lines.  If  the  ion  drag  and  viscosity  were  the  only  forcing  mechanisms 
involved,  strong  zonal  winds  should  prevail  once  the  plasma  is  lifted  along  the  field 
lines.  Likewise,  if  the  plasma  is  forced  down  the  field  lines,  the  zonal  winds  should  be 
reduced  because  of  the  increase  of  the  plasma  density.  Further  review  of  Figures  18  and 
19  reveals  an  inconsistency  in  the  December  solstice  meridional  and  zonal  winds. 
Before  midnight,  the  meridional  winds  are  in  a  northerly  direction  that  should  result  in 
forcing  the  plasma  down  the  field  lines,  but  the  zonal  winds  are  increasing  during  this 
time.  Figures  18  and  19  are  not  consistent  with  the  theory  outlined  above,  but  the 
meridional  wind  effect  on  increasing  or  decreasing  the  ion  drag  is  an  oversimplification 
of  the  physics  in  this  region.  This  theory  neglects  tidal,  gravity  wave,  and 
electrodynamic  effects.  This  raises  a  major  weakness  with  the  TIGCM93  model,  i.e., 
the  lack  of  modeling  the  electrodynamic  effects.  During  the  December  solstice,  a 
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vertical  ExB  drift  may  contribute  to  raising  the  plasma  and  indirectly  cause  the  strong 
eastward  winds.  This  lack  of  electric  field  predictions  prompted  the  creation  of  an 
enhanced  global  circulation  model,  TIE-GCM,  that  includes  the  electric  field 
calculations  [Richmond  et  al.,  1992],  We  do  not  use  TIE-GCM  predictions  since  this 
model  is  still  in  the  process  of  model  validation. 
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CHAPTER  4 

MEASUREMENT  TECHNIQUES  AND  DATA 

4.1.  INTRODUCTION 

Our  study  uses  measured  Fabry-Perot  Interferometer  (FPI)  and  Incoherent 
Scatter  Radar  (ISR)  winds.  In  addition  to  the  wind  data,  we  use  ISR  measured  electron 
density  data  for  calculating  the  airglow  volume  emission  rate  and  emission  layer  height 
This  chapter  provides  a  brief  outline  of  airglow  emissions,  FPI  equipment  and  operation, 
FPI  data  analysis,  ISR  neutral  wind  derivation,  ISR  electron  density  measurements, 
volume  emission  rate  calculations,  and  sphne  fits. 

4.2.  AIRGLOW  EMISSIONS 

Because  the  FPI  is  essentially  a  spectroscopic  device,  this  instrument  detects  and 
measures  spectral  lines.  Earth's  atmosphere  contains  many  naturally  emitting  sources  at 
altitudes  ranging  from  70  km  to  800  km.  These  emitting  sources  are  called  airglow. 
The  Earth's  airglow  results  from  the  re-emission  of  energy  from  atmospheric  particles 
following  daytime  excitation  by  solar  radiation  [Bone,  1991].  During  the  day,  the 
airglow  intensity  is  much  smaller  than  the  intensity  of  the  sunlight,  hence  making  it 
almost  impossible  to  detect.  At  night,  the  airglow  intensity  is  too  low  to  monitor  with 
the  unaided  eye,  but  this  low  intensity  presents  no  problem  for  a  spectroscopic 
instrument  such  as  the  FPI.  The  next  three  sections  highlight  the  most  commonly  used 
airglow  emissions  in  the  upper  atmosphere. 
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4.2.1.  Hydroxyl  (OH) 

The  OH  airglow  emission  occurs  over  a  relatively  narrow  region  of  the 
mesosphere.  The  peak  emission  occurs  around  85  km  and  the  thickness  of  this  emitting 
layer  is  approximately  10  km.  The  following  chemical  reaction  is  responsible  for  the 
OH  emission: 

O3  +  H  OH*  +  O2  (22) 

The  *  indicates  that  the  OH  molecule  is  in  the  excited  state  which  subsequently  emits  a 
photon  as  it  transitions  to  the  ground  state.  Although  OH  airglow  bands  cover  a  wide 
range  of  the  electromagnetic  spectrum,  the  infrared  OH  bands  are  the  strongest  emitters. 
The  OH  emission  intensity  ranges  from  values  greater  than  100  R  at  night  to 
approximately  10^  R  during  the  day  (where  1  R  =  10®  photons  cm'^sec'*).  The  OH 
emission  layer  is  important  since  it  characterizes  the  dynamics  in  the  mesosphere. 


4.2.2.  Monatomic  Oxygen,  0(^S),  557.7  nm 

At  night,  airglow  0(*S)  is  produced  in  both  the  lower  E  (97  km)  and  F  (230  km) 
ionospheric  regions.  The  dominant  processes  that  create  the  0(^S)  emissions  are: 

O  +  O  +  O-^OjH-  0(^S)  (Chapman  reaction)  (23) 


and 


0  +  0  +  M-^02+M 
o;  +  O  O2  +  0(‘S) 


(Barth  reactions) 


(24) 


The  *  identifies  molecules  in  the  excited  state  and  M  stands  for  any  third  body. 
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Typical  night  intensities  for  the  0(‘S)  airglow  range  from  100  R  to  200  R.  This 
airglow  has  a  radiative  lifetime  of  0.7  sec  [Jacka,  1984],  The  0(*S)  emission  line 
provides  a  method  to  study  the  E  and  F  ionospheric  region  dynamics. 

4.2.3.  Monatomic  Oxygen ,  0(^D),  630.0  nm 

Although  the  0(‘S)  airglow  emission  occurs  at  altimdes  near  250  km,  a  more 
commonly  used  emission  for  this  altitude  is  the  0(*D)  630.0  nm  airglow  emission.  The 
dissociative  recombination  of  O2  and  electrons  produce  the  0(^D)  airglow  emission. 
This  recombination  occurs  near  250  km  and  is  described  chemically  as: 

O2  +  c  — ^  0+0  (25) 

The  0(^D)'s  radiative  lifetime  of  110  sec  is  much  longer  than  the  0.7  sec 
radiative  lifetime  for  the  0(*S)  emission.  The  long  radiative  lifetime  and  long  colhsion 
time  between  F-region  particles  (0.3  sec)  make  the  0(*D)  emission  a  suitable  choice  for 
wind  and  kinetic  temperature  measurements  for  the  upper  thermosphere.  For  the 
remainder  of  this  thesis,  the  0(*D)  emission  line  tvill  be  referred  to  as  the  630.0  nm 
airglow  emission. 


4.3.  FABRY-PEROT  INTERFEROMETER 

4.3.1.  Historical  Perspective 

The  first  FPI  was  developed  in  1897  by  Charles  Fabry  and  Alfred  Perot, 
although  the  physical  principles  behind  the  FPI  were  actually  developed  more  than  a 
century  prior  by  Sir  Isaac  Newton  [Hernandez,  1986],  The  original  intent,  as  outlined 
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by  Fabry  and  Perot,  was  to  use  the  FPI  to  provide  accurate  measurements  of  the 
standard  meter.  It  was  not  long,  however,  before  interested  parties  realized  the 
spectroscopic  capabilities  of  this  instrument 

Astronomers  were  the  first  scientific  group  to  take  advantage  of  the  FPI's  high 
spectral  resolution.  High  spectral  resolution  is  necessary  to  determine  the  composition 
of  distant  stars.  Solar  physicists  used  the  instrument  to  characterize  the  sun's  spectrum 
and  discovered  other  elements  existed  in  the  solar  atmosphere  in  addition  to  hydrogen 
and  helium. 

In  1912,  scientists  used  the  FPI  to  test  the  kinetic  theory  of  gases.  Such  tests 
were  instrumental  in  explaining  the  Doppler  broadening  and  shift  effects  [Hernandez, 
1986].  Doppler  broadening  is  the  temperature  broadening  of  a  spectral  line  emission 
from  a  gas  compound  and  Doppler  shift  is  the  wavelength  shift  as  a  result  of  the  bulk 
motion  of  the  emitters.  Both  Doppler  broadening  and  doppler  shift  are  essential  when 
describing  the  upper  atmosphere.  Additional  uses  of  the  FPI  include  measurements  of 
viscosity  effects  in  air,  the  Zeeman  effect,  the  index  of  refraction  of  gases,  and  the 
isotropic  and  fine  structure  investigations  [Hernandez,  1986]. 

The  first  upper  atmospheric  wind  and  temperature  FPI  measurements  were 
conducted  in  1960.  Since  then,  the  FPI  has  routinely  been  used  to  measure  the  winds 
and  temperatures  of  the  mesosphere  and  thermosphere. 

4.3.2.  Interferometry 

This  section  outlines  the  basic  principles  of  interferometry.  Optics  textbooks 
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by  Hecht  and  Zajac  [1979]  and  Fowles  [1975]  present  a  more  complete  description  of 
this  subject. 

4.3.3.  Two-beam  and  Multiple-beam  Interference 

AU  interferometic  instruments,  including  the  FPI,  are  based  on  the  principle  of 
optical  interference.  A  thorough  description  of  optical  interference  requires  quantum 
mechanical  explanations,  an  endeavor  beyond  the  scope  of  this  study.  Fortunately,  the 
theory  behind  the  FPI  requires  only  basic  optical  principles. 

To  understand  interference,  it  is  helpful  to  picture  a  light  wave  as  a  sinusoidal 
wave.  Each  light  wave  has  a  unique  amplitude  (intensity),  wavelength,  and  relative 
phase  angle.  For  simplicity,  we  consider  light  of  only  one  wavelength  (monochromatic 
light). 

If  two  light  waves  are  directed  at  the  same  point  in  space  but  travel  along 
different  paths,  the  waves  wiU  affect  each  other's  amplitude  either  constructively  or 
destructively.  Because  we  imposed  a  monochromatic  condition  on  these  light  waves, 
the  resultant  amplitude  from  the  interference  of  these  two  beams  is  a  function  of  the 
relative  phase  angle  difference  between  the  two  light  waves.  This  type  of  interference  is 
called  two-beam  interference. 

Two-beam  interference  occurs  for  coherent  and  incoherent  light  Incoherent 
light  consists  of  light  waves  with  different  amplitudes  and  relative  phase  angles.  An 
example  of  incoherent  light  is  the  scattered  light  from  the  earth's  atmosphere.  This  light 
undergoes  a  near  continuous  interference  process  but  it  is  impossible  to  observe  due  to 
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the  different  phase  angles  and  the  high  frequencies  (10^^  Hz).  Any  light  interference  in 
earth's  atmosphere  occurs  too  quickly  to  be  detected  by  human  eyes.  Coherent  light  is 
light  whose  waves  have  a  constant  phase  difference.  It  is  easy  to  observe  interference 
with  coherent  light  because  it  is  possible  to  find  a  distance  where  light  interferes  and  the 
resultant  wave  amplitude  remains  constant.  It  is  coherent  light  that  produces  the  bright 
(constructive)  and  dark  (destructive)  fringes  seen  in  most  interference  patterns. 

Multiple-beam  interference  differs  from  two-beam  interference  in  the  number  of 
light  waves  (or  light  beams)  that  interfere.  The  beams  in  multiple-beam  interference  are 
made  up  of  more  than  two  light  waves  with  constant  phase  differences  and  identical 
wave  amplitudes.  This  type  of  interference  occurs  for  the  transmitted  beams  from  two 
parallel,  semireflective  mirrors.  These  transmitted  beams  are  parallel  to  each  other  and 
theoretically  interfere  at  infinity.  To  observe  this  interference,  a  positive  lens  is  placed 
in  the  path  of  the  parallel  beams,  thus  causing  the  light  beams  to  constructively  or 
destructively  interfere  at  the  lens's  focal  length.  It  is  multiple-beam  interference  that 
forms  the  optical  foundation  for  the  FPL 

4.3.4.  Basic  Optical  Principles  for  the  FPI 

The  most  common  method  of  producing  a  larger  number  of  mutually  coherent 
beams  is  by  division  of  wave  amplitude  [Fowles,  1975].  Two  parallel,  semitransparent 
mirrors  are  aU  that  is  needed  to  produce  this  effect.  A  light  beam  with  amplitude  Eq 
incident  on  the  outside  mirror  surface  will  be  reflected  between  the  two  mirrors  plates 
(see  Figure  20).  Figure  20  shows  the  resulting  parallel  transmitted  waves  (Eot ,  Eot  r  , 
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etc.,  where  t  and  r  are  the  transmission  and  reflection  coefficients,  respectively). 

The  phase  difference  between  two  successive  transmitted  waves  is  given  by: 

5  =  ^  nd  cosQ  (26) 

where  8  is  the  phase  difference  between  the  two  rays,  Xq  is  the  wavelength  of  the 
incident  radiation,  n  is  the  index  of  refraction  of  the  medium,  d  is  the  mirror  separation, 
and  0  is  the  angle  of  incidence.  Figure  21  shows  a  detailed  look  at  the  successive  waves 
and  the  geometric  factors  used  in  the  phase  difference  calculation.  From  Figure  21,  the 
path  difference  between  two  successive  transmitted  waves,  t,  is: 

t  =  2ndcosQ  (27) 

The  light  waves  in  that  emerge  from  the  second  mirror  are  parallel  to  each  other. 
This  situation  is  similar  to  the  transmitted  beams  discussed  in  the  previous  section.  If  a 


Eoth^ 


E.t^ 


Figure  20.  Paths  of  light  rays  in  multiple  reflection  between  two  parallel  mirrors  [After 
G.  Fowels,  1975,  pg.  86.  Reproduced  with  permission  of  Dover  Publications] 
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positive  lens  is  placed  beyond  the  second  mirror,  the  rays  will  be  brought  together  at  the 
focal  plane  of  the  lens.  At  this  point,  enough  information  has  been  provided  to 
understand  the  most  basic  FPI  arrangement  (see  Figure  22).  Figure  22  shows  how  the 
focusing  lens  brings  the  rays  together  at  a  specific  point.  The  etalon  in  the  diagram  is 
identical  to  the  two  parallel  mirror  setup  discussed  earlier  in  this  study. 


Figure  21.  Diagram  showing  the  path  difference  between  two  successive  rays  [After  G. 
Fowles,  1975,  pg.  87.  Reproduced  with  permission  of  Dover  Publications] 


Broad  Collimating  Focusing  Focal  Ring  system 

source  lens  lens  plane 


Figure  22.  Basic  Fabry-Perot  arrangement  [After  G.  Fowles,  1975,  pg.  91. 
Reproduced  with  permission  of  Dover  Publications] 
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The  ring  system  in  Figure  22  is  produced  when  a  broad  source  of  Hght  enters  the 
FPL  Since  airglow  is  a  broad  light  source,  the  incident  light  enters  the  FPI  over  a  wide 
range  of  angles  [Vadnais,  1992],  After  the  light  is  transmitted  by  the  etalon,  the 
focusing  lens  focuses  the  light  at  the  focal  plane  where  the  ring  system  is  formed.  Each 
ring  in  the  ring  system  represents  a  different  order  of  interference.  The  bright  rings  are 
created  by  constructive  interference  and  the  dark  rings  by  destructive  interference. 

At  the  angle  ([)  where  a  bright  ring  or  fringe  is  located,  the  relationship  between 
the  wavelength  and  the  mirror  (etalon)  spacing  d  is: 

mX  =  2nd  cos  ^  (28) 

where  m  is  an  integer  identifying  the  order  of  interference  [Gomall,  1983].  If  the  light 
rays  interfere  on  the  optical  axis  ((])  =  0),  equation  (28)  reduces  to: 

mX  =  2nd  (29) 

An  FPI  satisfying  equation  (29)  may  be  tuned  to  transmit  any  wavelength  by  vaiying  the 
index  of  refraction  n  or  the  mirror  spacing  d.  The  index  of  refraction  may  be  changed 
by  varying  the  gas  pressure  between  the  mirrors.  This  pressure  scanning  method  was 
used  to  collect  the  Arequipa  and  Arecibo  FPI  data  used  for  this  study.  Today,  most 
FPIs  scan  by  changing  the  mirror  spacing  d.  The  most  versatile  technique  for  changing 
the  mirror  spacing  is  using  piezoelectric  elements.  These  elements  keep  the  mirrors 
parallel  to  one  another  as  the  spacing  changes  [Gomall,  1983]. 

The  intensity  of  the  transmitted  waves  is  an  important  quantity  for  FPI 
measurements.  The  equation  for  the  intensity  of  the  transmitted  rays  is: 
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(30) 

where  A  is  the  absorptance,  R  the  reflectivity,  5  the  phase  difference,  and  F  the  finesse. 
The  term 

1  +  Fsin"-  (31) 

2 

is  known  as  the  Airy  function.  The  finesse  F  is  discussed  in  the  next  section. 


4.3.4. 1.  The  Meaning  of  Finesse 

Finesse  is  a  measure  of  the  interferometer's  abihty  to  resolve  closely  spaced 
lines;  the  higher  the  finesse  the  better  [Gomall,  1983].  A  more  conceptual  definition  of 
finesse  is  the  number  of  interfering  monochromatic  hght  beams  involved  in  forming  the 
interference  fringes.  Many  factors  such  as  the  reflectivity,  mirror  flatness,  and  pinhole 
diameter  affect  the  total  instrument  finesse  value. 

The  reflectivity  flnesse  for  a  plane  mirror  interferometer  is  a  function  of  the 
mirror  reflectivities  [Gomall,  1983].  The  equation  for  reflectivity  finesse  is: 


(32) 


where  i?  is  the  reflectivity. 

No  mirror  is  perfectly  flat  A  non-flat  mirror  used  in  an  FPI  becomes  important 
when  viewing  small  spectral  resolution.  The  equation  for  this  flatness  flnesse  is: 


(33) 
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where  M  is  defined  as  the  fractional  wavelength  deviation  from  true  flatness  or 
parallelism  across  the  mirror  aperture  [Gomall,  1983]. 

The  net  finesse  due  to  flatness  and  reflectivity  is  called  the  instrumental  finesse 
and  is  mathematically  defined  as: 


(34) 


When  the  FPI  is  used  as  a  spectrometer,  the  pinhole  size  determines  the  degree 
of  coUimation  of  hght  passing  and  reaching  the  detector  [Gomall,  1983].  If  the  pmhole 
size  is  too  large,  many  hght  rays  pass  through  the  system.  The  effect  of  too  many  rays 
entering  the  FPI  detector  is  to  broaden  the  instrumental  line  width.  This  pinhole  finesse 
is  defined  as: 


D^d 


(35) 


where/  is  the  focal  length  of  the  second  lens  and  £)  is  the  pinhole  diameter. 

The  total  instrumental  finesse  for  a  FPI  spectrometer  is  a  function  of  the 
reflectivity,  flatness,  and  pinhole  finesse.  Mathematically  the  total  finesse  for  the  FPI 
spectrometer  is  defined  as: 


(36) 


4.3.4.2.  Throughput  and  Etendue 

The  strength  of  the  FPI  over  other  interferometers  is  its  "etendue"  or  light 
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gathering  power  [Gomall,  1983].  For  the  ideal  case  of  small  aperture  and  flat  parallel 
mirrors,  the  transmission  on  the  peak  of  a  fringe  is: 


where  A  is  the  absorption  due  to  the  mirrors.  This  maximum  transmission  results  when 
the  phase  difference,  5,  is  a  multiple  of  2ji.  The  etendue  for  the  FPI  is  mathematically 
defined  as: 

U  =  Cl.A^=:!^  (38) 

^  4dF. 

where  X  is  the  incident  wavelength,  O  is  the  solid  angle  subtended  at  a  mirror  aperture 
Af^,  D  is  the  pinhole  diameter,  d  is  the  distance  between  the  mirrors,  and  Ff  is  the 
instrument  finesse. 

4.3.4.3.  Etalon 

A  glass  plate  with  optically  flat  and  parallel  surfaces  separated  by  a  fixed 
distance  is  the  simplest  form  of  Fabry-Perot  and  is  called  an  etalon  [Higgins,  1987].  A 
light  ray  incident  to  the  first  parallel  surface  will  undergo  many  reflections  between  the 
two  parallel  mirror  surfaces.  If  the  light  is  coherent,  constructive  and  destructive 
interference  will  occur  between  the  reflected  and  transmitted  portions  of  the  beam 
[Higgins,  1987]. 

4.3.5.  FPI  Equipment  Description 

This  section  provides  a  basic  understanding  of  theArequipa  and  Arecibo  FPIs. 
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We  do  not  intend  to  provide  a  detailed  description  of  the  FPI  equipment  and  operation, 
but  we  want  to  provide  enough  information  to  understand  how  the  FPI  measures  the 
upper  atmospheric  winds. 

High-resolution  determinations  of  the  center  frequency  of  the  oxygen  line 
emission  of  630.0  nm  with  reference  to  the  center  frequency  in  zenith  observations 
(averaged  over  a  night  of  vertical  measurements)  give  a  direct  indication  of  the 
component  of  thermospheric  flow  in  the  observing  direction  [Meriwether  et  al,  1986]. 

4.3.5. 1.  ArequipaFPI 

A  detailed  description  of  the  Arequipa  FPI  is  provided  in  Biondi  et  al.  [1990  and 
1991]  and  Meriwether  et  al.  [1986].  The  Arequipa  measurements  were  made  at  the 
Arequipa  Automated  Airglow  Observatory  (16.5°S,  71.5°W,  3.2°S  dip  latitude)  in  Peru. 
The  Arequipa  observatory  is  located  ~  450  km  south  of  Jicamarca,  Peru. 

A  refractive-index-tuned,  field-widened  FPI  of  94-mm  clear  aperture  and 
2x10”^  nm  resolution  is  used  to  determine  the  Doppler  shifts  in  the  nightglow  line 
emissions  over  Arequipa.  The  complete  FPI  system  consists  of  the  optical  instrument 
and  a  companion  tilting  filter  photometer. 

The  wind  measurements  are  taken  along  the  cardinal  azimuths  at  30°  elevation 
and  in  the  zenith,  with  the  average  630.0  nm  line  position  looking  vertically  used  to 
define  the  zero  velocity  reference  for  the  eastward  and  northward  wind  component 
determination.  Measurements  of  the  632.8  nm  line  from  a  stabilized  He-Ne  laser  are 
used  to  determine  the  free  spectral  range  and  overall  finesse  of  the  FPI,  as  well  as  any 
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instrument  drifts  [Biondi  et  al,  1990].  The  spectral  profile  of  the  source  emission  was 
scanned  by  changing  the  index  of  refraction  of  the  medium  (argon  gas)  between  the 
plates  by  compressing  or  expanding  tlie  argon  with  a  computer-controlled  piston  housed 
in  a  sealed  cylinder  [Meriwether  etal.,  1986]. 

The  observations  and  calculations  indicate  that  the  typical  airglow  emission  is 
broad  (>  50  km  between  the  half-intensity  points)  and  that  the  altitude  of  peak  emission 
varies  between  ~  250  km  and  ~  325  km,  depending  of  the  time  of  night,  season,  and  part 
of  the  solar  cycle.  It  is  assumed  that  the  630.0  nm  winds  do  not  change  appreciably  in 
magnitude  and  direction  over  this  altitude  range. 

4.3.5.2.  AreciboFPI 

The  Arecibo  FPI  is  a  pressure-scanned  FPI  with  a  15-cm  diameter  etalon 
situated  near  the  town  of  Arecibo,  Puerto  Rico  (18°N,  67°W,  50°  dip  latitude).  The 
pressure  of  the  argon  gas  in  the  system  is  controlled  by  a  piston  that  is  stepped  to 
various  positions  under  computer  control.  The  field  of  view  of  the  FPI  is  1/2°.  The 
aperture  size  is  6.3  mm,  a  size  which  reduces  the  overall  instrument  finesse  to  about  10. 
A  50-mm  diameter  interference  filter  with  a  bandwidth  of  .5  nm  was  used  to  restrict 
transmission  to  some  25  orders  near  630.0  nm.  The  etalon  spacing  was  10  mm, 
minimizing  the  contamination  of  OH  emissions. 

At  Arecibo,  FPI  measurements  are  made  at  eight  equally  spaced  azimuth 
positions  with  a  constant  zenith  angle  of  60°.  In  analyzing  the  data,  it  is  assumed  that 
the  spatial  variation  of  the  neutral  wind  velocity  is  linear  across  the  region  that  is 
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observed.  The  eight  independent  measurements  allow  calculations  of  the  deformation 
of  the  flow,  the  mean  velocity,  and  the  divergence. 

Because  the  mean  overhead  neutral  wind  velocity  is  determined  by  observing  the 
airglow  from  eight  emitting  regions,  each  located  some  500  km  from  the  station, 
systematic  errors  are  possible  unless  a  number  of  conditions  hold.  These  conditions  are 
the  horizontal  gradient  in  wind  velocity  is  constant  across  the  observing  region,  the 
height  of  the  emitting  layer  is  the  same  everywhere,  and  the  neutral  wind  does  not  vary 
with  altitude. 

4.3.6.  FPI  Measured  Winds 

The  next  two  sections  present  the  Arequipa  and  Arecibo  data  measurements. 
4.3.6. 1.  Arequipa 

TABLES  4  and  5  show  a  breakdown  of  the  Arequipa  eastward  measured  FPI 
data  in  terms  of  10.7-cm  solar  flux  (O)  and  magnetic  activity  (given  by  the  Kp  index). 
The  total  data  points  and  the  percentage  of  the  total  are  given  for  each  25-unit  solar  flux 
increment  TABLE  4  shows  that  the  1988-1989  data  are  characterized  by  moderate  to 
high  solar  activity.  Although  the  northward  data  are  not  shown,  these  data  show  similar 
total  and  percentage  of  the  total  values. 

TABLE  5  shows  the  total  eastward  measured  FPI  data  points  for  a  given  Kp 
interval  and  a  given  3-hour  period.  This  table  shows  the  majority  of  data  were  collected 
during  quiet  and  moderate  magnetic  activity  (Kp  < 

Figure  23  shows  the  15-day  (1/2  month)  distribution  for  the  measured  Arequipa 
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TABLE  4.  Solar  Flux  Histogram  of  1988-1989  Arequipa  Neutral  Wind  Data 
Total  Number  of  Wind  Profiles:  474 


Part  (a) 


Solar  Flux 

<  100 

100  <  O  <  125 

125  <  O  <  150 

Total 

0 

50 

68 

Percentage 

of 

Total  (%) 

0.0 

10.5 

14 .  i 

Part  (b) 

Solar  Flux 

150  < ^><  175 

175  <  O  <  200 

<I)>200 

Total 

106 

120 

130 

Percentage 

of 

22.4 

25.3 

27.4 

Total  (%) 

TABLE  5.  Magnetic  Activity  (Kp)  Histogram  of  1988-1989  Arequipa  Neutral  Wind 

Data 


Part(a) 


Hour  (LT) 

0  <  Kp  <  1 

1  <  Kp  <  2 

2  <  Kp  <  3 

3  <  Kp  <  4 

<  21 

68 

201 

181 

91 

21-00 

153 

497 

619 

310 

00-03 

165 

452 

688 

264 

03-06 

65 

226 

241 

76 

65 


data.  This  graph  shows  the  majority  of  Arequipa  data  were  collected  during  the  May- 
August  months;  hence,  the  winter  solstice  is  weU  represented.  The  least  amount  of  data 
was  collected  during  the  November-February  months  (summer  solstice).  The 
northward  15-  day  distribution  is  not  shown  but  is  similar  to  the  results  shown  in  Figure 
23. 

4.3.6.2.  Arecibo 

TABLES  6  and  7  show  a  breakdown  of  the  Arecibo  eastward  measured  FPI 
data  in  terms  of  10.7-cm  solar  flux  and  magnetic  activity.  TABLE  6  shows  the  Arecibo 
data  are  represented  for  low,  moderate,  and  high  solar  flux  activity.  TABLE  7  shows 
the  majority  of  the  data  were  collected  for  quiet  magnetic  activity  (Kp  <  3’). 

The  15-day  distribution  for  the  measured  Arecibo  data  is  presented  in  Figure  24. 
This  graph  shows  the  majority  of  Arecibo  data  were  collected  during  May-AugusL 
Unlike  the  Arequipa  data,  the  Arecibo  data  show  a  wide  monthly  distribution  allowing  a 
more  complete  seasonal  average  study.  The  northward  15-day  Arecibo  distribution  is 
not  shown  but  is  similar  to  the  results  shown  in  Figure  24. 

4.3.7.  FPI  Data  Analysis 

Seasonal  averages  of  the  nightly  variations  in  the  eastward  and  northward  winds 
have  been  obtained  from  the  measurements  by  taking  aU  the  data  points  for  a  given 
season  for  a  selected  10.7-cm  solar  flux  and  magnetic  activity  intervals,  sorting  them 
into  30-minute  bins,  and  averaging  the  number  of  data  points,  n,  in  each  bin  to 
obtain  the  sample  mean  value,  =(5^Vi)/n.  The  variability  in  the  measurements 


Number  of  Days  with  Data 
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AREQUIPA 

Fabry- Perot  15-Day  Distribution  (looking  East  &  West) 


Month  (15- Day  Increments) 


Figure  23.  Frequency  distribution  of  Arequipa  data  collection.  The  time  interval  is  15- 
day  (1/2  month)  increments. 


TABLE  6.  Solar  Flux  Histogram  of  1981-1989  Arecibo  Neutral  Wind  Data 
Total  Number  of  Wind  Profiles:  202 


Part  (a) 


Solar  Flux 

<  100 

100  <  O  <  125 

125  <  O  <  150 

Total 

64 

17 

29 

Percentage 

of 

Total  (%) 

31.7 

8.4 

14.4 

Part(b) 


Solar  Flux 

150  < <I)<  175 

175  <  d)<  200 

O>200 

Total 

33 

32 

27 

Percentage 

of 

Total  (%) 

16.3 

15.8 

13.4 

TABLE  7.  Magnetic  Activity  (Kp)  Histogram  of  1981-1989  Arecibo  Neutral  Wind 

Data 


Part  (a) 


Hour  (LT) 

0  <  Kp  <  1 

1  <  Kp  <  2 

2  <  Kp  <  3 

3  <  Kp  <  4 

<  21 

50 

104 

148 

71 

21-00 

139 

250 

365 

195 

00-03 

146 

368 

327 

159 

03-06 

137 

284 

259 

70 

Part  (b) 


Hour  (LT) 

4  <  Kp  <  5 

5  <  Kp  <  6 

6  <  Kp  <  7 

Kp  > 

<  21 

32 

7 

3 

2 

72 

17 

6 

7 

00-03 

92 

17 

9 

8 

55 

34 

4 

0 
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ARECIBO 


Month  (15-Day  Increments) 


Figure  24.  Frequency  distribution  of  Arecibo  data  collection.  The  time  interval  is  15- 
day  (1/2  month)  increments. 
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for  each  30-minutB  bin  is  determined  by  the  sample  mean  standard  deviation,  s,  of  the 
data  from  its  mean.  The  quantity  s  is  given  by 

For  our  study,  we  analyze  eastward  and  northward  wind  using  graphs  like  those 
in  Figures  25  and  26.  These  plots  are  separated  into  three  parts:  the  bottom  plot  shows 
the  scatter  profiles  of  each  night’s  wind  observations,  the  top  histogram  indicates  the 
number  of  points  in  each  30-minute  time  bin,  and  the  middle  are  the  sample  mean  wind 
component  values  at  30-minute  intervals  (open  circles)  and  the  ±1  standard  deviations 
from  the  means  (error  bars).  The  middle  plot  also  shows  the  10.7-cm  solar  flux  and  the 
magnetic  activity  intervals,  the  representative  day  for  the  season,  and  the  mean  solar  flux 
(4>)  and  magnetic  activity  (Kp )  values. 

Figure  25,  the  May- August  eastward  wind  over  Arequipa,  shows  large  standard 
deviations  (~  ±  75  m/s  -  100  m/s).  We  use  a  statistical  analysis  technique  outlined  in 
Biondi  et  al  [1991]  to  show  that  we  can  apply  a  reasonable  confidence  limits  in  spite  of 
the  large  nightly  wind  variation.  Figure  27  shows  the  northward  and  eastward  data  for 
the  given  time,  solar  flux,  and  magnetic  activity  conditions  follow  the  form  of  a  normal 
distribution.  Because  the  data  follow  a  normal  distribution,  we  can  apply  the  usual 
confidence  limit  analysis  to  provide  us  a  more  quantitative  assessment  in  the  average 
wind  profile.  TABLE  8  shows  the  results  of  apply  the  95%  confidence  limits  to  the  data 
in  Figure  25.  TABLE  8  shows  the  small  lower  and  upper  boundaries  of  the  confidence 
interval,  which  provides  us  more  confidence  in  the  average  wind  profile  for  Arequipa. 


Eastward  Wind  (m/s)  Eastward  Wind  (m/s)  number  of  points 
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Arequipa  Fabry- Perot  1988,1989,1990 
Averages  for  May- August  looking  East  +  West 
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Figure  25.  May-August  eastward  average  of  the  thermospheric  wind  over  Arequipa. 
The  top  plot  is  the  number  of  data  points  in  each  30  min  interval,  the  middle  plot  is  the 
mean  values  of  the  wind  component  (points)  and  the  ±  1  standard  deviation  (vertical 
bars),  and  the  bottom  is  the  scatter  plots  of  the  eastward  wind  for  the  given  solar  flux 
and  magnetic  activity  intervals. 
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Northward  Wind  (m/s)  Northward  Wind  (m/s)  number  of  points 
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Arecibo  Fabry- Perot 

Averages  for  November- February  looking  North  +  South 


Universal  Time 

Figure  26.  November-February  northward  average  of  the  thermospheric  wind  over 
Arecibo.  The  top  plot  is  the  number  of  data  points  in  each  30  min  interval,  the  middle 
plot  is  the  mean  values  of  the  wind  component  (points)  and  the  ±  1  standard  deviation 
(vertical  bars),  and  the  bottom  is  the  scatter  plots  of  the  eastward  wind  for  the  given 
solar  flux  and  magnetic  activity  intervals. 


Number  of  Eastward  Wind  Points  Number  of  Northward  Wind  Points 


'2 


Arequipa  Fabry-Perot  (1988,1989,1990) 
Averages  for  May-August 


Figure  27.  May- August  frequency  distribution  bar  graphs  of  the  deviation  of  the  wind 
values  in  each  bin  from  their  mean  (v-Vmean),  for  Arequipa.  The  smooth  curves  are  the 
least-squares  fits  to  the  data  of  normal  distributions  centered  at  (v-Vmcan)  =  0. 
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TABLE  8.  95%  Confidence  Intervals  for  the  May- August  Eastward  Winds  over 

Arequipa,  Peru. 


Time 

(UT) 

Vmean 

(m/s) 

Standard 

Deviation 

(m/s) 

Number  of 
Points 

95% 

confidence 

lower 

boundary 

(m/s) 

95% 

confidence 

upper 

boundary 

(m/s) 

-0.75 

172 . 1 

125.8 

60 

140.2 

203.9 

-0.25 

162.1 

106.7 

97 

140.8 

183.3 

0.25 

162.8 

108.3 

75 

138.3 

187.3 

0.75 

159.1 

80.7 

76 

141.0 

177.3 

1.25 

175.5 

88.7 

98 

157.9 

193.0 

1.75 

184.4 

84.3 

135 

170.2 

198.6 

2.25 

183.1 

82.7 

172 

170.8 

195.5 

2.75 

185.4 

79.7 

202 

174.4 

196.4 

3.25 

181.4 

75.5 

219 

171.4 

191.4 

3.75 

176.1 

84.1 

228 

165.2 

187.1 

4.25 

164.8 

76.2 

232 

155.0 

174.6 

4.75 

157.8 

77.6 

217 

147.5 

168.1 

5.25 

139.9 

74.1 

219 

130.1 

149.7 

5.75 

138.6 

76.9 

197 

127.8 

149.3 

6.25 

129.5 

76.5 

184 

118.4 

140.5 

6.75 

116.2 

75.8 

167 

104.7 

127.7 

7.25 

111.8 

80.1 

145 

98.7 

124.8 

7.75 

106.2 

85.7 

151 

92.6 

119.9 

8.25 

99.8 

90.9 

115 

83.2 

116.4 

8.75 

102.0 

89.7 

91 

83.6 

120.4 

9.25 

104.7 

80.8 

69 

85.7 

123.8 

9.75 

90.9 

93.0 

73 

69.5 

112.2 

10.25 

23.1 

99.2 

75 

0.6 

45.5 
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We  did  not  conduct  a  similar  analysis  for  Arecibo  since  the  Arecibo  standard  deviations 
about  the  mean  are  small  compared  to  those  of  Arequipa. 

4.4.  ISR  NEUTRAL  WIND  MEASUREMENTS 

In  addition  to  using  FPI  measured  winds,  we  use  neutral  winds  derived  from  the 
ISR  at  Arecibo.  These  ISR  winds  provide  us  additional  data  for  comparing  with  the  FPI 
measured  meridional  winds  over  Arecibo.  We  can  compare  directly  the  ISR-derived 
neutral  winds  to  the  meridional  FPI  measured  winds  since  the  magnetic  declination  at 
Arecibo  is  less  than  10°. 

4.4.1.  ISR  Neutral  Winds 

The  determination  of  the  meridional  neutral  wind  velocity  along  the  geomagnetic 
meridian  from  ISR  data  is  a  well  understood  technique  [e.g.,  Harper,  1973;  Salah  and 
Holt,  1974;  Burnside  et  al,  1983,  1991;  Wickwar  et  al,  1984].  For  our  study,  we  use 
the  Arecibo  ISR-derived  neutral  winds  from  1983  to  1987.  A  complete  description  of 
the  ISR  neutral  wind  derivation  at  Arecibo  is  outlined  in  Burnside  et  al.  [1983].  In  this 
section,  we  present  a  brief  description  of  the  steps  used  to  calculated  the  neutral  wind 
from  the  Arecibo  ISR  data. 

Data  are  obtained  by  the  Arecibo  430  MHz  ISR  in  a  continuous  beamswinging 
mode  with  a  fixed  zenith  angle  of  15°.  Calculations  for  the  F-layer  parameters  are 
obtained  every  16  minutes  (the  time  it  takes  to  complete  one  rotation  of  the  feed). 
Line-of-sight  plasma  velocities  are  derived  from  the  imaginary  part  of  the 
autocorrelation  function.  Next,  a  harmonic  analysis  is  performed  on  the  line-of-sight  ion 
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velocities  to  get  the  horizontal  components.  For  these  calculations,  it  is  assumed  that 
spatial  and  temporal  gradients  in  the  ion  velocity  are  small  enough  to  be  neglected. 

In  the  direction  of  geomagnetic  south,  the  neutral  wind  velocity  is  given  by 
Burnside  et  al.  [1983] 

=  (v,  -  Vj )  sec  7  (40) 

where  vn  is  the  ion  velocity  antiparallel  to  the  magnetic  field,  Vd  is  the  O'^  diffusion 


velocity,  and  7  is  the  magnetic  inclination  angle  (~  50°).  The  diffusion  velocity  of  in 
the  neutral  atmosphere  was  derived  by  Schunk  and  Walker  [1970]  and  is  given  by 


Vj  =  -D„  sin  7 


1  I  dN  1  dT  w  dT 


-+■ 


+- 


l_77p  N  dz  Tp  dz  dz 


(41) 


where  N  is  the  electron  density,  =  IkT^lmg  is  the  plasma  scale  height,  Tp  =  (Ti  + 
Te)  /  2,  Tr  =  (Ti  +  Tn)  /  2,  T„  is  the  neutral  temperature,  and  w  is  a  thermal  diffusion 
factor  with  a  value  of  0.37.  The  parameter  m  is  the  ion  mass  (16  amu),  and  g  is  the 
acceleration  due  to  gravity.  The  ambipolar  diffusion  coefficient.  Da,  is  given  by  Schunk 
and  Walker  [1970] 


3.02x10^^ 

a[0]  +  l9.9[N^]T;^^^ 


(42) 


where 

a  =  1.08-0.141og,o  r„  +0.0045(logio  T„f  (43) 

The  ISR  measurements  of  vn  are  obtained  every  38  km  and  the  determination  of 
Vd  require  knowledge  of  the  plasma  scale  height  and  the  vertical  gradient  in  electron 
density,  which  is  obtained  from  the  600-m  height  resolution  Barker  code  measurements 
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at  Arecibo  [Burnside  et  al,  1991].  Except  for  the  ambipolar  diffusion  coefficient.  Da, 
each  of  the  five  terms  in  equation  (41)  for  Vd  is  obtained  by  ISR  measurements. 
Measurements  of  atomic  oxygen  [O]  and  nitrogen  [N2]  are  found  from  the  MSIS  model. 
These  concentrations  are  used  for  calculating  Dgin  equation  (42). 

4.4.2.  ISR  Measured  Data 

The  Arecibo  ISR  data  used  in  this  research  were  obtained  from  the  National 
Center  for  Atmospheric  Research  (NCAR)  Couphng,  Energetics,  and  Dynamics  of  the 
Atmospheric  Regions  (CEDAR)  database.  The  CEDAR  database  is  located  at  the 
NCAR  facility  in  Boulder,  Colorado.  For  our  study,  we  use  collected  data  from  1983  to 
1987;  hence,  this  data  characterize  the  neutral  wind  during  low  solar  activity.  TABLES 
9  and  10  show  the  breakdown  of  this  ISR  data  in  terms  of  solar  activity  and  magnetic 
activity.  Figure  28  shows  an  example  of  the  ISR  neutral  wind  graphs  used  for  our 
study. 


4.5.  ELECTRON  DENSITY  DATA 

We  use  Jicamarca,  Peru  and  Arecibo  measured  ISR  electron  density  profiles  for 
calculating  the  airglow  emission  heights  over  Arequipa  and  Arecibo,  respectively. 
Jicamarca  is  450  km  north  of  Arequipa;  hence  the  measurements  over  Jicamarca  are  also 
representative  of  the  conditions  over  Arequipa.  The  electron  density  profiles  were 
obtained  from  the  CEDAR  database.  TABLES  11  and  12  show  the  breakdown  of  the 
solar  flux  and  magnetic  activity  for  Jicamarca  and  TABLES  13  and  14  show  the 
breakdown  for  the  Arecibo  electron  density  data. 
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TABLE  9.  Solar  Flux  Histogram  of  1983-1987  ISR  Neutral  Wind  Data 
Total  Number  of  Data  Points;  100 

Part  (a) _ 


Solar  Flux 

>  100 

100  < ^><  125 

125  <  O  <  150 

Total 

59 

22 

19 

Percentage 

of 

Total  (%) 

59 

22 

19 

Part  (b) 


Solar  Flux 

150  <  Ox  175 

175  <  O  <  200 

^>>200 

Total 

0 

0 

0 

Percentage 

of 

Total  (%) 

0 

0 

0 

TABLE  10.  Magnetic  Activity  (Kp)  Histogram  of  1983-1987  ISR  Neutral  Wind  Data 
Part  (a) 


Part  (b) 


Meridional  Wind  (m/s)  Meridional  Wind  (m/s)  number  of  points 
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Arecibo  ISR  Winds  1983,1984,1985,1986,1987 
Averages  for  November- February 


20  1  6  11  16 


20  1  6  11  16 
Universal  Time 


Figure  28.  November-February  northward  average  of  the  ISR  calculated  neutral  wind 
over  Arecibo.  The  top  plot  is  the  number  of  data  points  in  each  1  hr  interval,  the  middle 
plot  is  the  mean  values  of  the  wind  component  (points)  and  the  ±  1  standard  deviation 
(vertical  bars),  and  the  bottom  is  the  scatter  plots  of  the  eastward  wind  for  the  given 
solar  flux,  magnetic  activity,  and  altitude  intervals. 
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TABLE  11.  Solar  Flux  Histogram  of  1966-1969, 1984, 1988,  and  1991  Jicamarca 

Electron  Density 


Part  (a) 


Solar  Flux 

>100 

100  < <1><  125 

125  <  O  <  150 

Percentage 

of 

Total  (%) 

1.5 

8.3 

52.8 

Part  (b) 

Solar  Flux 

150  <  ^  <  175 

175  <  O  <  200 

0>2OO 

Percentage 

of 

Total  (%) 

23.1 

8.7 

5.6 

TABLE  12.  Magnetic  Activity  (Kp)  Histogram  of  Jicamarca  Electron  Density  Data 

Percentage  of  Total 


Part  (a) 


Hour  (LT) 

0  <  Kp  <  1 

1  <  Kp  <  2 

2  <  Kp  <  3 

3  <  Kp  <  4 

00-04 

20.3 

26.0 

23.9 

20.1 

04-08 

18.0 

24.0 

27.0 

18.6 

08-12 

18.1 

28.4 

18.8 

17.3 

12-16 

15.6 

29.7 

25.3 

11.1 

16-20 

15.9 

28.5 

26.6 

12.7 

20-00 

16.3 

24.8 

25.4 

21.1 

Part  (b) 


Hour 

(LT) 

4  <  Kp  <  5 

5  <  Kp  <  6 

6  <  Kp  <  7 

Kp  >  7 

00- 

-04 

8.5 

1.0 

0.2 

0.0 

04- 

-08 

8.4 

2.5 

1.1 

0.4 

08- 

-12 

8.8 

4.9 

2.2 

1.4 

12- 

-16 

11.9 

2.7 

1.6 

2.0 

16- 

-20 

10.0 

3.6 

1.5 

1.2 

20- 

-00 

8.6 

1.9 

1.1 

0.7 
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TABLE  13.  Solar  Flux  Histogram  of  1981-1989  Arecibo  Electron  Density 


Part  (a) 


Solar  Flux 

>100 

100  < «I><  125 

125  <  <I>  <  150 

Percentage 

of 

40.3 

12.5 

14.7 

Total  (%) 

Part  (b) 


Solar  Flux 

150  <  <f»  <  175 

175  <  ^)<  200 

0>200 

Percentage 

of 

11.2 

6.8 

14.4 

Total  (%) 

TABLE  14.  Magnetic  Activity  (Kp)  Histogram  of  Arecibo  Electron  Density  Data 

Percentage  of  Total 


Part  (a) 


Hour  (LT) 

0  <  Kp  <  1 

1  <  Kp  <  2 

2  <  Kp  <  3 

3  <  Kp  <  4 

00-04 

9.1 

27.7 

29.2 

19.7 

04-08 

12.1 

28.7 

26.9 

22.2 

08-12 

12.0 

24.6 

32.7 

19.0 

12-16 

11.7 

25.9 

30.5 

17.9 

16-20 

12.3 

19.0 

31.6 

20.4 

20-00 

10.8 

19.0 

28.8 

23.1 
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4.5.1.  Airglow  Emissions  and  Heights 

The  ionosphere-thermosphere  is  home  to  several  important  chemical  reactions. 
For  this  study,  we  focus  our  attention  on  the  dissociative  recombination  of  Oj.  This 
reaction  plays  a  fundamental  role  in  the  generation  of  two  of  the  most  important 
emission  lines  widely  used  for  FPI  studies,  that  is,  the  557.7  nm  and  630.0  nm  emissions 
[Sobral  et  al,  1993].  The  dissociative  recombination  of  Oj  is  the  main  source  of  the 
atomic  oxygen  excited  species  0(^D)  and  0(^S)  in  the  nocturnal  F-region  [Sobral  et  al, 
1993].  The  equations  and  formulas  used  to  calculate  the  volume  emission  rates  are 
found  in  Link  and  Cogger  [1988]. 


4.5.2.  Calculating  the  Volume  Emission  Rate 

TABLE  15  lists  the  chemical  reactions  and  their  associated  rate  coefficients  used 
in  determining  the  630.0  nm  emission  rate.  The  rate  coefficient  or  rate  of  reaction  is 
defined  as  the  decrease  in  concentration  of  a  reactant  or  the  increase  in  concentration  of 
a  product  in  a  unit  time.  The  electron,  ion,  and  neutral  temperatures  are  given  by  the 
symbols  Te,  Ti,  and  r„,  respectively.  Last,  =  7^  1300  and  X-,=TJ  300. 

The  630.0  nm  volume  emission  rate  is  given  by  the  following  formula: 


y  0.76Pifc.[On[O.] 

“  l+(A:3[iV2]4-^,[02]-Hfc5[e])/A„ 

where  Pj  =  1.1  and  =  7.45  x  lO'^s”^. 


(44) 


The  Jicamarca  and  Arecibo  electron  density  profiles  have  a  38-km  height 
resolution.  To  compensate  for  the  large  38-kra  resolution,  we  used  a  numerical 
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TABLE  15.  Volume  Emission  Rate  Parameters 

Reaction  Rate  Coefficient,  (cm' ^s'^)  References 

O"  +  O2  O;  +  0  =  Chen  etal.  [1978] 

0('D)  +  N2  -4  0('P)  +  N2  =2.0xl0~"e‘“  ®^”  Streitetal.  [1976] 

OOD)  +  O2  ->  0(^P)  +  O2  ;t4=2.9xl0""e^-^^”  Streitetal.  [1976] 

0('D)  +  e  0(^P)  +  e  ^'5  =1.6x10"’^  7^°^^  Berrington  and  Burke 

_ [1981] _ 

technique  to  fit  the  known  electron  density  data  points  to  a  5*  degree  polynomial.  The 
numerical  method,  least-squares  fitting  of  discrete  points,  is  described  in  detail  in 
Hombeck  [1983].  After  obtaining  the  5"*  degree  polynomial,  we  impose  a  minimum 
altitude  boundary  of  170  km,  a  maximum  altitude  boundary  of  530  km,  and  a  5  km 
altitude  increment.  Figure  29  shows  the  polynomial  fit  (dashed  line)  and  the  average 
(solid  line)  electron  density  profile.  Once  the  electron  density  profile  is  determined  by 
the  polynomial,  we  calculate  the  volume  emission  rate. 

As  indicated  in  equation  (44),  we  must  determine  the  concentration  of  O2, 
and  N2  before  calculating  the  630.0-nm  volume  emission  rate.  We  obtain  the  0“^  from 
the  electron  density  profile  and  the  neutral  concentrations  from  the  MSISE90  model. 
The  details  and  description  of  this  model  are  outlined  in  Hedin  [1991].  The  required 
inputs  for  the  MSISE90  model  are  the  solar  flux  and  the  magnetic  activity  given  by  the 
Ap  index.  For  these  inputs,  we  enter  the  average  solar  flux  and  the  Ap  index  provided  of 
the  electron  density  data  for  a  given  solar  activity  and  magnetic  activity  interval. 

The  time  interval  for  the  data  analysis  is  +  15  minutes  either  side  of  a  given  hour 
(e.g.,  we  use  the  data  within  ±  15  minutes  of  0100  LT  to  determine  the  average 
electron  density  profile  for  0100  LT). 
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ARECI  BO  81 , 82, 83, 84, 85, 86, 87, 88, 89 

Averages  for  May- August  Hour:  0200UT 


630.  Onm  Volume  Emission  Rate  (photons  cm'^sec'^) 


Figure  29.  May-August  calculated  630.0  nm  volume  emission  rate  over  Arecibo.  The 
top  plot  is  calculated  volume  emission  rate,  the  middle  plots  are  the  mean  electron 
density  profile  (solid  line)  and  the  modeled  electron  density  profile  (dashed  line),  and  the 
bottom  is  the  scatter  plots  of  the  eastward  wind  for  the  given  solar  flux  and  magnetic 
activity  intervals. 
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4.5.3  F-Layer  Height  and  Airglow  Emission  Height 

Using  the  electron  density  profiles  measured  at  Arecibo  and  Jicamarca, 
MSISE90  empirical  model  predictions,  and  rate  coefficients  and  formulas  defined  in 
Link  and  Cogger  [1988  and  1989],  we  calculate  the  volume  emission  rates.  These 
volume  emission  rates  wiU  provide  us  the  airglow  emission  height  or  the  height  at  which 
the  FPI  measures  the  thermospheric  winds.  Since  we  are  interested  in  establishing  the 
climatological  season  average,  we  calculate  the  airglow  emission  heights  for  all  seasons, 
high  and  low  solar  activity,  and  quiet  and  disturbed  magnetic  activity. 

4.5.3. 1.  Arecibo  F-Layer  and  Airglow  Emission  Heights 

As  described  in  section  4.5.3,  we  use  the  Arecibo  electron  density  profiles 
collected  during  January  1981  to  November  1989  for  our  study.  This  data  set  covers 
both  high  and  low  solar  activity  and  quiet  and  disturbed  magnetic  activity  effects. 
Figures  30  and  26  show  the  solar  cycle  and  magnetic  activity  effects  on  the  F-layer  peak 
(Hmax)  and  airglow  emission  height  (h^MAx)  for  all  seasons. 

The  solar  activity  effects  (see  Figure  30)  were  separated  for  high  (^>  >  120)  and 
low  (O  <  120)  conditions  for  quiet  magnetic  activity  (Kp  <  3)  for  the  winter  solstice 
(November-February),  equinox  (March-April,  September-October),  and  summer 
solstice  (May- August). 

For  the  winter  solstice,  we  note  h,,MAX  is  ~  260  km  and  235  km  through  the  night 
for  high  and  low  solar  activity,  respectively.  The  Hmax  is  ~  75  km  and  85  km  higher 
than  hrjMAx  for  the  pre-  and  postmidnight  periods  for  high  flux  conditions  while 
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Hmax  is  ~  75  km  higher  than  h^MAx  through  the  night  for  low  flux  conditions. 

The  equinox  months  for  high  solar  activity  show  hriMAx  at  ~  275  km  in  the 
premidnight  period  lowering  to  255  km  in  the  postmidnight  period  and  hT,MAx  is  ~  250 
km  and  240  km  in  the  pre-  and  postmidnight  periods  for  low  solar  activity.  The  Hmax  is 
~  120  km  and  100  km  higher  than  ^max  for  the  pre-  and  postmidnight  periods  for  high 
flux  conditions  and  Hmax  is  ~  75  km  higher  than  h,,MAx  for  low  solar  activity. 

The  summer  solstice  months  show  h,,MAx  at  ~  280  km  in  the  premidnight  period 
and  ~  260  km  in  the  postmidnight  period  for  high  solar  activity  and  ~  250  km  and  ~  240 
km  in  the  pre-  and  postmidnight  periods.  Hmax  is  ~  120  km  and  100  km  higher  than 
hi^MAx  for  the  pre-  and  postmidnight  periods  during  high  solar  activity  while  Hmax  is  ~ 
80  km  higher  than  h^MAX  for  the  night  at  low  flux  conditions. 

The  magnetic  activity  effects  (see  Figure  31)  were  separated  for  disturbed  (Kp  > 
3)  and  quiet  (Kp  <  3)  conditions  for  all  solar  flux  activity.  We  note  no  significant 
difference  between  the  quiet  and  disturbed  hriMAx  for  aU  seasons.  We  further  note  that 
the  average  ^max  height  is  250  km.  Hmax  also  shows  no  significant  difference  between 
the  quiet  and  disturbed  conditions  for  all  seasons.  Hmax  shows  the  largest  difference  in 
the  equinox  and  winter  months  (~  120  km  higher  than  ^max  in  the  premidnight  period) 
and  smallest  difference  in  the  winter  solstice  months  (~  75  km  and  100  km  higher  than 
hijMAx  in  the  pre-  and  postmidnight  periods) 

4.5.3.2.  Jicamarca  F-Layer  and  Airglow  Emission  Heights 

The  Jicamarca  electron  density  data  were  collected  in  1966-1969,  1984, 1988, 


ARECIBO  F- Layer  and  Airglow  Emission  Profile 
A  -  H .  n  -  h .  All  Fluxes 


Nov-  Feb 


=  4.0 
Kp  =  1.6 


400 


lar- Apr,  Sep- Oct 

- =  3.8 

^ 

Kp  =  1.8 

May-Aug 

- =  3.8 

Kp  =  2.0 

Local  Time 

Figure  31.  Arecibo  seasonal  F-layer  and  airglow  emission  height  averages  for  quiet  and 
disturbed  magnetic  activity. 
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and  1991.  These  data  were  characterized  by  low  to  moderate  solar  activity.  Figures  32 
and  33  show  the  solar  cycle  and  magnetic  activity  effects  on  the  Hmax  and  hT,MAx  for  aU 
seasons. 

Figure  33  shows  the  moderate  (d>  >  130)  and  low  (<&  <  130)  solar  activity 
effects  on  Hmax  and  ^max  for  low  magnetic  activity  conditions  for  the  summer  solstice 
(November-February),  equinox  (March-April,  September-October),  and  winter  solstice 
(May-August).  Due  to  the  limited  Jicamarca  data  set,  it  was  necessary  to  increase  the 
magnetic  activity  interval  to  Kp  <  4  to  ensure  a  good  statistics  for  the  moderate  and  low 
solar  activity  average. 

For  the  summer  solstice,  we  note  h,,MAx  is  ~  230  km  and  250  km  in  the  pre-  and 
postmidnight  periods  for  the  moderate  flux  activity  while  hi^MAX  is  ~  270  km  in  the 
premidnight  period  for  low  solar  activity.  The  Jicamarca  data  set  lacked  sufficient  data 
satisfying  the  criteria  outline  above  for  an  Hmax  and  h^MAx  in  the  postmidnight  period 
for  low  flux  activity.  The  Hmax  is  ~  300  km  higher  than  ^max  in  the  premidnight  period 
and  100  km  higher  than  h^MAx  in  the  postmidnight  period  for  moderate  flux  conditions 
while  Hmax  is  ~  150  km  higher  than  h^MAx  in  the  premidnight  period  for  low  solar 
activity.  The  large  Hmax  values  in  the  premidnight  period  is  a  result  of  the  prereversal 
enhancement  of  the  ExB  drift  after  sunset  over  the  magnetic  equator  [Anderson, 
1974]. 

The  equinox  months  for  high  solar  activity  show  h,,MAx  at  ~  225  km  in  the 
premidnight  period  lowering  to  240  km  in  the  postmidnight  period  for  moderate  and 
low  solar  activity.  The  Hmax  is  ~  300  km  and  90  km  higher  than  hr^MAx  for  the  pre-  and 
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postmidnight  periods  for  moderate  and  low  solar  flux  condition. 

The  winter  solstice  months  show  ^max  at  ~  250  km  through  the  night  for 
moderate  and  low  solar  activity.  Hmax  is  ~  175  km  and  75  km  higher  than  ^max  for  the 
pre-  and  postmidnight  periods  during  moderate  and  low  solar  activity. 

The  magnetic  activity  effects  (see  Figure  33)  were  separated  for  disturbed  (Kp  > 
3)  and  quiet  (Kp  <  3)  conditions  for  all  solar  flux  conditions.  Similar  to  the  results  for 
Arecibo,  the  data  show  no  significant  variation  between  the  quiet  and  disturbed  hi^MAx 
for  all  seasons.  For  the  summer  and  winter  solstice  months,  we  note  an  average  hi^MAx 
of  ~  250  km  while  the  equinox  months  show  an  hi^MAx  of  -220  km  and  -  240  km  in  the 
pre-  and  postmidnight  periods.  Hmax  also  shows  no  significant  difference  between  the 
quiet  and  disturbed  conditions  for  the  equinox  and  winter  solstice  months,  but  we  do 
note  a  lower  Hmax-  (~  50  km)  for  quiet  magnetic  activity  during  the  summer  solstice 
months.  Last,  aU  seasons  show  a  significantly  larger  Hmax  when  compared  to  ^max  for 
the  premidnight  period  for  the  summer  and  equinox  months. 

4.6.  SPLINES 

We  use  an  empirical  data  modeling  technique  using  spline  fits  outlined  in  Wand 
[1981]  to  model  the  F-region  neutral  winds  over  Arequipa  and  Arecibo.  Splines  are 
bell-shaped  functions  whose  degree  of  continuity  depends  on  the  order  of  the  spline. 
We  use  splines  of  degree  four,  which  means  these  splines  have  continuous  derivatives 
up  to  the  second  order.  We  use  this  spline  technique  because  it  provides  good 
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agreement  with  the  average  wind  profile,  thereby  making  the  spline  result  an  excellent 
choice  for  comparison  with  modeled  data. 

We  constructed  two  separate  spline  models  of  the  thermospheric  wind  data,  one 
for  the  zonal  component  of  the  neutral  wind  and  the  other  for  the  meridional 
component.  Nodes  for  these  models  were  set  every  1.5  hours  for  a  total  of  nine  nodes. 
We  found  this  node  distribution  enhanced  the  smoothness  of  the  spline  fit  curves  to  the 
measured  results. 

The  zonal  and  meridional  components  of  the  neutral  wind  models  are  both 
defined  by  the  form; 


1 


a,- 1  +a,_2  sm 


f2nd\  (2nd\ 

-  +a,-,cos - 

l,365j  U65j 


+ai^flux + +  Ui^flux^ 


N,{t) 


(45) 


where  m  is  the  number  of  nodes  in  the  respective  model,  and  aij  are  the  model 
coefficients  obtained  by  a  least  squares  fit  to  the  respective  data  points.  Because  of  the 
differences  between  the  eastward  and  northward  wind  speeds,  the  final  empirical  model 
will  be  different  for  the  two  wind  components  because  the  coefficients,  fly,  will  be 
different  The  model  given  in  equation  (45)  is  a  function  of  the  local  time  it),  the  10.7- 
cm  solar  flux  (flux),  the  magnetic  activity  (K^,  and  the  day  of  the  year  (d).  We  perform 
least  squares  fits  to  the  respective  zonal  and  meridional  data  sets  to  obtain  the  model 
coefficients.  Figures  34  to  36  show  examples  of  the  agreement  between  the  spline  fits 
and  the  average  wind  profiles. 


Arequipa  Fabry- Perot  1988,1989,1990 
Averages  for  May- August  looking  East 


West 


O  172.9 

Day  180 

1  <  0<  400 

Kp  2.2 

0.0  <  Kp  <  9.9 

Universal  Time 

Figure  34.  Spline  fit  for  the  May-August  eastward  wind  over  Arequipa.  The  solid  line 
in  the  middle  plot  is  the  spline  fit  for  the  average  solar  flux  and  magnetic  activity 
conditions. 


Northward  Wind  (m/s)  Northward  Wind  (m/s)  number  of  points 


Arecibo  Fabry- Perot 

Averages  for  November- February  looking  North  +  South 


Universal  Time 

Figure  35.  Spline  fit  for  the  Noveraber-February  northward  wind  over  Arecibo.  The 
solid  line  in  Ae  middle  plot  is  the  spline  fit  for  the  average  solar  flux  and  magnetic 
activity  conditions. 


Meridional  Wind  (m/s)  Meridional  Wind  (m/s)  number  of  points 
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Arecibo  I  SR  Winds  1983,1984,1985,1986,1987 
Averages  for  November- February 


20  1  6  11  16 
Universal  Time 


Figure  36.  Spline  fit  for  the  November-February  Arecibo  ISR  calculated  neutral  wind 
along  the  magnetic  meridian.  The  solid  line  in  the  middle  plot  is  the  spline  fit  for  the 
average  solar  flux  and  magnetic  activity  conditions. 
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CHAPTER  5 

QUIET  TIME  LOW  LATITUDE  THERMOSPHERIC  WINDS 

5.1.  INTRODUCTION 

This  chapter  presents  the  geomagnetic  quiet  measured  wind  data  over  Arequipa, 
Peru  and  Arecibo,  Puerto  Rico.  We  calculated  the  average  winds  during  quiet  time 
conditions  to  eliminate  magnetic  activity  effects.  In  addition  to  Fabry-Perot  wind 
measurements,  we  use  ISR  calculated  neutral  winds  over  Arecibo  to  enhance  our  study 
of  seasonal  and  solar  cycle  effects  on  the  neutral  winds.  We  also  investigate  the 
TIGCM93  and  HWM93  predicted  winds  during  quiet  geomagnetic  conditions. 

5.2.  SEASONAL  AND  SOLAR  CYCLE  EFFECTS 

This  section  details  the  seasonal  and  solar  cycle  effects  on  the  thermospheric 
winds  over  Arequipa  and  Arecibo.  With  the  measured  data,  we  also  use  TIGCM93  and 
HWM93  results  to  provide  an  indication  of  the  success  of  these  models  in  predicting  the 
quiet  time  average  winds  over  these  low  latitude  locations. 

5.2.1.  Arequipa 

To  analyze  both  the  seasonal  and  solar  cycle  effects  on  the  neutral  winds,  we  use 
1988-1990  measured  FPI  data.  In  all  the  following  Arequipa  FPI  data  analysis,  we 
neglect  the  December  solstice  (November-February)  measurements  due  to  the  limited 
data  collected  during  this  period. 
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5.2.1. 1.  Seasonal  Effects 

Figures  37  and  38  show  the  averaged  eastward  and  northward  F-region  quiet 
time  (Kp  <  3)  thermospheric  neutral  winds  over  Arequipa  for  all  solar  flux  levels  during 
equinox  (March-April,  September-October)  and  June  solstice  (May-August)  periods. 
These  average  zonal  winds  in  Figure  37  are  generally  eastward  through  most  of  the 
night  for  aU  seasons  except  after  ~  0300  LT  where  the  data  show  an  eastward  decrease 
or  westward  reversal.  The  strongest  winds  occur  during  the  June  solstice  months  with 
average  peak  winds  of  175  to  250  m/s  in  the  preraidnight  period. 

Figure  38  shows  the  average  meridional  winds  over  Arequipa.  These  average 
winds  consist  of  small  wind  speeds  (<  50  m/s)  after  2100  LT  for  the  equinox  and  June 
solstice  while  the  winds  before  2100  LT  show  a  significant  northward  wind  speed  (>  50 
m/s)  for  the  equinox  months  and  a  southward  wind  (>  50  m/s)  for  the  June  solstice 
period. 

Despite  using  a  different  analysis  technique,  our  results  are  consistent  with 
previous  Arequipa  wind  studies  outlined  in  Biondi  et  al.  [1991].  Biondi  et  al.  [1991] 
investigated  the  seasonal  effects  by  calculating  monthly  and  seasonal  wind  averages 
(equinox  and  June  solstice)  for  a  given  year  (1988,  1989,  etc.).  Our  averaging  method 
differs  from  that  used  in  Biondi  et  al.  [1991]  because  we  combined  the  1988-1990  wind 
data  for  a  given  season.  We  also  used  measured  FPI  winds  during  quiet  geomagnetic 
conditions  to  reduce  magnetic  activity  effects,  whereas  Biondi  et  al.  [1991]  combined 
quiet  and  disturbed  winds  for  the  seasonal  average  calculations.  Biondi  et  al.  [1991] 
showed  the  strongest  eastward  winds  occur  during  the  June  solstice  period  and  the 
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Figure  37.  Measured  FPI  eastward  wind  averages  over  Arequipa,  Peru,  for  equinox 
(March-April,  September-October)  and  near  June  solstice  (May-August)  for  all  solar 
activity  and  quiet  (Kp  <  3)  geomagnetic  activity  conditions. 
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Figure  38.  Measured  FPI  northward  wind  averages  over  Arequipa,  Peru,  for  equinox 
and  near  June  solstice  for  all  solar  activity  and  quiet  geomagnetic  activity  conditions. 
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northward  winds  exhibit  moderate  (<  50  m/s)  wind  speeds  after  midnight  during  the 
equinox  and  June  solstice  periods. 

We  also  compared  our  equinox  meridional  wind  results  to  derived  neutral  winds 
from  the  Atmospheric  Explorer  E  (AE-E)  satellite.  AE-E  measurements  of  the  630.0 
run  emission  in  the  nighttime  equatorial  thermosphere  were  used  to  infer  the  height  of 
the  F-region  peak,  Hmax,  during  the  December  solstice  and  spring  equinox  (March- 
April)  periods  [Burrage  et  al,  1990].  These  airglow  measurements  were  collected 
during  February  1980  through  April  1981,  a  period  characterized  by  moderate  to  high 
solar  flux  activity  and  low  to  moderate  geomagnetic  activity  [Burrage  et  al,  1990]. 
Using  the  calculated  Hmax  values,  Burrage  et  al  [1990]  obtained  magnetic  meridional 
wind  components  as  a  function  of  dip  latitude  and  local  time.  Our  equinox  meridional 
wind  average  does  not  agree  with  the  AE-E-derived  wind  average  for  the  spring 
equinox  period  at  the  lowest  latitudes.  The  AE-E-derived  wind  near  the  magnetic 
equator  shows  a  southward  wind  reversal  beginning  near  2100  LT  until  0230  LT  while 
our  FPI  result  for  equinox  shows  a  northward  wind  from  2100  LT  to  0430  LT.  Further 
review  of  the  AE-E  winds  shows  the  derived  winds  from  15°S  to  25°S  dip  latitudes  are 
consistent  with  our  equinox  seasonal  average  wind.  Burrage  et  al  [1990]  stated  that 
the  derived  winds  at  the  lowest  latitudes  are  subject  to  the  greatest  uncertainties  since 
the  contribution  of  the  meridional  neutral  wind  to  the  change  in  Hmax  decreases  to  zero 
at  the  equator. 

5.2. 1 .2.  Solar  Cycle  Effects 

Figures  39  and  40  show  the  averaged  eastward  and  northward  F-region  quiet 
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time  thermospheric  neutral  winds  over  Arequipa  for  two  levels  of  the  10.7-cm  solar  flux 
(O)  index  during  equinox  and  June  solstice  periods.  The  1988-1990  Arequipa  FPI  data 
are  characterized  by  moderate  (d>  <  160)  to  high  («I>  >  160)  solar  flux  activity.  The 
open  symbols  represent  the  mean  wind  speed  for  the  given  solar  activity  and  the  error 
bars  represent  +  1  standard  deviation  about  this  mean.  The  error  bars  for  the  high  solar 
flux  means,  the  closed  symbols,  are  not  shown  but  they  are  comparable  to  those  of  the 
moderate  solar  flux  conditions. 

Figure  39  shows  mainly  eastward  winds  for  all  seasons  over  Arequipa.  The  data 
show  no  significant  variation  due  to  solar  activity  during  the  equinox  period;  however, 
the  data  show  a  systematic  eastward  increase  in  the  postmidnight  period  from  moderate 
to  high  solar  flux  conditions  during  the  June  solstice  months.  Figure  40  shows  no 
significant  variation  of  the  neutral  meridional  wind  due  to  solar  cycle  effects. 

Biondi  et  al  [1990  and  1991]  reported  that  the  nighttime  zonal  wind  increases 
with  solar  activity.  These  previous  studies  showed  larger  eastwards  in  the  premidnight 
sector  from  solar  minimum  to  maximum  while  our  results  show  the  larger  eastward 
winds  in  the  postmidnight  period  for  similar  solar  activity  conditions.  One  possible 
explanation  of  the  discrepancy  between  our  results  and  the  previous  studies  is  that 
Biondi  et  al.  [1990  and  1991]  used  FPI  measurements  during  1983-1986  and  1988-1990 
while  we  limited  our  data  set  to  1988-1990.  The  1983-1986  FPI  data  are  characterized 
by  low  solar  flux  conditions.  These  data  also  showed  lower  eastward  winds  in  the 
premidnight  period  compared  to  the  1988-1990  zonal  wind  data  set  We  did  not  use  the 
1983-1986  wind  data  set  in  our  analysis  because  the  FPI  instrument  encountered  several 
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Figure  39.  Measured  eastward  wind  averages  over  Arequipa,  Peru,  for  various  seasons 
for  moderate  and  high  solar  activity.  The  average  neutral  wind  velocities  are  obtained 
by  examining  1/2  hour  time  bins  for  the  given  solar  flux  interval  (160  <  4>  <  160)  and 
quiet  geomagnetic  conditions.  The  triangle  symbol  represents  bin  averages  consisting  of 
less  than  10  data  points. 
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Figure  40.  Measured  northward  wind  averages  over  Arequipa,  Peru,  for  various 
seasons  for  moderate  and  high  solar  activity.  The  average  neutral  wind  velocities  are 
obtained  by  examining  1/2  hour  time  bins  for  the  given  solar  flux  interval  (160  <  <^  < 
160)  and  quiet  geomagnetic  conditions.  The  triangle  symbol  represents  bin  averages 
consisting  of  less  than  10  data  points. 
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problems  during  these  years.  Our  northward  wind  results  agree  with  the  previous 
meridional  wind  studies  in  Biondi  et  al.  [1991]. 

5.2. 1.3.  TIGCM93  Solar  Cycle  Activity  Model  Results 

Figures  41  and  42  show  the  TIGCM93  solar  cycle  activity  predictions  for 
Arequipa  for  all  seasons.  To  properly  investigate  the  solar  cycle  effects,  we  analyze  the 
wind  predictions  for  high  (<I>  =  195),  moderate  (<I>  =  135),  and  low  (<I>  =  75)  solar 
activity  conditions.  We  use  the  predicted  quiet  time  winds,  Kp  =  1,  to  reduce  magnetic 
activity  effects. 

Figure  41  shows  the  eastward  predictions  for  the  December  solstice  (day  =  360), 
equinox  (day  =  90),  and  June  solstice  (day  =  180).  Similar  to  the  measured  FPI  winds 
in  Figure  39,  the  modeled  winds  are  generally  eastward  throughout  the  night  for  all 
seasons  with  a  westward  reversal  occurring  after  0430  LT  for  days  360  and  90.  The 
TIGCM93  peak  eastward  winds  are  generally  50  -  75  m/s  less  than  the  measured  FPI 
average  winds.  Exact  agreement  between  the  predicted  and  measured  wind  does  not 
determine  the  success  of  the  TIGCM  prediction.  The  relative  difference  between  the 
predicted  low  and  high  solar  activity  winds  compared  to  the  relative  difference  for  the 
measured  winds  provides  a  better  measure  of  how  the  TIGCM93  predicts  the 
experimental  wind  results. 

We  compare  only  the  average  zonal  winds  for  the  equinox  and  June  solstice 
periods  to  the  TIGCM  predictions.  The  measured  or  experimental  data  set  does  not 
contain  enough  moderate  solar  flux  data  during  November  -  February  for  a  valid 
comparison  during  the  December  solstice  period.  Because  the  measured  Arequipa  data 
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Figure  41.  TIGCM93  eastward  wind  predictions  over  Arequipa,  Peru,  for  various 
seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet  magnetic 
activity. 
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Figure  42.  TIGCM93  northward  wind  predictions  over  Arequipa,  Peru,  for  various 
seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet  magnetic 
activity. 
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are  characterized  by  moderate  to  high  solar  activity,  we  use  the  predicted  high  (O  = 
195)  and  moderate  (<I>  =  135)  relative  wind  difference  to  compare  to  the  measured 
data.  For  day  90,  the  TIGCM93  predicts  a  westward  velocity  perturbation  from  high  to 
moderate  flux  conditions  in  the  premidnight  period  while  the  postmidnight  shows  no 
significant  variation  due  to  solar  activity.  The  relative  differences  between  the  predicted 
wind  results  for  the  premidnight  period  are  not  consistent  with  the  experimental  results 
while  the  postmidnight  model  results  show  better  agreement  with  the  experimental 
moderate  and  high  flux  wind  difference.  For  day  180,  the  model  results  show  no 
significant  wind  variation  due  to  solar  activity.  The  T1GCM93  predictions  agree  with 
the  experimental  wind  results  in  the  premidnight  period  but  show  less  agreement  in  the 
postmidnight  period.  During  the  postmidnight  period,  the  measured  winds  show  a 
systematic  eastward  increase  from  moderate  to  high  solar  activity  while  the  model 
predicts  no  significant  variation  of  the  winds  due  to  solar  activity.  Figure  42  shows  no 
significant  variation  of  the  meridional  averages  due  to  solar  activity  effects,  a  result 
consistent  with  the  measured  winds. 

5.2. 1.4.  HWM93  Solar  Cycle  Activity  Model  Predictions 

Figures  43  and  44  show  HWM93  solar  cycle  activity  predictions  for  Arequipa 
for  all  seasons.  Similar  to  the  previous  section,  we  analyze  the  wind  predictions  for  high 
(O  =  195),  moderate  (<[>  =  135),  and  low  (O  =  75)  solar  activity  conditions.  We  use  the 
HWM93  predicted  quiet  time  winds,  Kp  =  2,  to  compare  to  the  FPI  measured  winds. 
This  Kp  value  agrees  with  the  average  Kp  index  for  the  measured  FPI  wind  under  similar 
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solar  flux  and  magnetic  activity  conditions. 

Because  the  HWM93  model  uses  FPI  measured  winds  for  Arequipa,  we  do  not 
expect  to  find  significant  variations  between  the  measured  and  predicted  results.  In 
addition  to  the  FPI  data,  the  HWM93  also  contains  ISR  and  satellite  data,  allowing  us 
to  obtain  model  predictions  during  times  of  limited  or  no  FPI  wind  measurements. 
Figure  43  shows  the  v^requipa  eastward  predictions  for  the  December  solstice  (Day  = 
360),  equinox  (Day  =  90),  and  June  solstice  (Day  =  180).  These  predicted  zonal  results 
for  the  equinox  and  June  solstices  are  in  complete  agreement  with  our  measured  data 
results  (see  Figure  39).  The  day  360  predicted  results  show  httle  variation  with  solar 
activity  in  the  premidnight  period  and  a  westward  perturbation  from  high  to  low  solar 
activity  in  the  postmidnight  period.  Figure  44  shows  the  model  predictions  for  the 
Arequipa  meridional  winds  are  not  affected  by  solar  cycle  variations;  a  result  that  agrees 
with  our  experimental  data. 

5.2.2.  Arecibo 

We  use  1980-1989  measured  FPI  data  to  analyze  both  the  seasonal  and  solar 
cycle  effects.  In  addition  to  the  measured  FPI  data,  we  use  ISR-derived  neutral  winds 
to  supplement  the  FPI  data  analysis.  In  the  following  Arequipa  FPI  solar  flux  analysis, 
we  neglect  the  equinox  data  analysis  because  of  the  limited  high  solar  flux  (<1>  >  130)  data 
collected  during  this  period. 

5.2.2. 1.  Seasonal  Effects 

Figures  45  and  46  show  the  averaged  eastward  and  northward  F-region  quiet 
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Figxme  43.  HWM93  eastward  wind  predictions  over  Arequipa,  Peru,  for  various 
seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet  magnetic 
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Figure  44.  HWM93  northward  wind  predictions  over  Arequipa,  Peru,  for  various 
seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet  magnetic 
activity. 
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time  thermospheric  neutral  winds  over  Arecibo  for  all  solar  flux  levels  during  December 
solstice,  equinox,  and  June  solstice  periods.  Figure  45  shows  that  the  averaged  FPI 
winds  are  generally  eastward  through  the  night  during  the  December  solstice.  The 
equinox  and  summer  solstice  show  westward  reversals  at  approximately  0300  LT  and 
2400  LT,  respectively.  The  maximum  eastward  winds  occur  during  the  December 
solstice  with  average  peak  wind  values  ranging  from  75  to  130  m/s  during  the 
premidnight  period. 

Figure  46  presents  the  average  meridional  winds.  These  winds  show  a 
pole ward-equatorward-pole ward  wind  flow  through  the  night  for  aU  seasons.  The 
maximum  meridional  winds  occur  during  the  June  solstice  where  the  data  show  peak 
southward  winds  of  75  to  100  m/s. 

Burnside  and  Tepley  [1989]  performed  similar  studies  of  the  seasonal  effects  of 
the  thermospheric  winds  over  Arecibo.  One  major  difference  between  this  previous 
study  and  our  work  is  Burnside  and  Tepley  [1989]  used  measured  FPI  data  from  1980- 
1987  while  our  study  uses  these  data  and  additional  FPI  measurements  during  1988  and 
1989.  Burnside  and  Tepley  [1989]  showed  that  the  strongest  eastward  winds  are 
observed  during  the  December  solstice  and  equinox  periods,  which  agrees  with  our  data 
set.  This  previous  study  also  shows  the  reversal  to  westward  flow  in  the  postmidnight 
period  during  the  June  solstice,  which  agrees  with  our  eastward  wind  results.  Our 
meridional  wind  averages  agree  with  the  results  outlined  in  Burnside  and  Tepley  [1989]. 
This  previous  study  showed  the  meridional  winds  are  weakest  during  the  December 
solstice  and  strongest  during  the  June  solstice  months.  Burnside  and  Tepley  [1989]  also 
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Figure  45.  Measured  eastward  wind  averages  over  Arecibo,  Puerto  Rico,  for  various 
seasons  for  all  solar  activity  and  quiet  magnetic  activity  conditions. 
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showed  the  strongest  and  weakest  meridional  wind  speeds  are  observed  during  the  June 
and  December  solstices,  respectively.  These  wind  directions  and  magnitudes  are 
consistent  with  our  meridional  results. 

5.2.22.  Solar  Cycle  Effects  (FPI  Measurements) 

Figures  47  and  48  show  the  averaged  eastward  and  northward  F-region  quiet 
time  thermospheric  neutral  winds  over  Arecibo  for  two  levels  of  the  10.7-cm  solar  flux 
index  for  the  December  solstice  and  June  solstice  periods.  These  data  were  separated 
for  low  solar  activity  (O  <  130)  and  high  solar  activity  (O  >  130). 

Figure  47  shows  generally  eastward  winds  throughout  the  night  for  the 
December  solstice.  The  averaged  zonal  wind  for  the  June  solstice  months  is  eastward  in 
the  premidnight  sector,  reverses  before  0100  LT,  and  becomes  westward  for  the  rest 
of  the  night.  These  data  show  a  small  westward  perturbation  for  most  of  the 
night  during  the  December  solstice  period  from  solar  minimum  to  maximum.  We  note 
no  significant  variation  between  the  low  and  high  solar  activity  averages  for  the  equinox 
and  June  solstice  periods.  Figure  48  shows  no  significant  variation  between  the  low  and 
high  solar  activity  northward  averaged  winds  during  the  December  solstice  and  June 
solstice  periods;  however,  the  data  show  perhaps  a  small  southward  velocity  perturbation 
for  the  equinox  period  from  solar  minimum  to  maximum. 

Burnside  and  Tepley  [1989]  concluded  that  the  thermospheric  neutral  winds 
over  Arecibo  are  unaffected  by  changes  in  the  solar  cycle.  Our  study  shows  solar  cycle 
effects  on  the  zonal  and  meridional  neutral  wind  averages.  We  enhanced  the  solar  cycle 
analysis  study  performed  in  Burnside  and  Tepley  [1989]  by  using  a  lower  time  bin 
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Figure  47.  Measured  eastward  wind  averages  over  Arecibo,  Puerto  Rico,  for  various 
seasons  for  low  and  high  solar  activity.  The  average  neutral  wind  velocities  are 
obtained  by  examining  1/2  hour  time  bins  for  the  given  solar  flux  interval  (130  <  <E>  < 
130)  and  quiet  geomagnetic  conditions.  The  triangle  symbol  represents  bin  averages 
consisting  of  less  than  10  data  points. 
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Figure  48.  Measured  northward  wind  averages  over  Arecibo,  Puerto  Rico,  for  various 
seasons  for  low  and  high  solar  activity.  The  average  neutral  wind  velocities  are 
obtained  by  examining  1/2  hour  time  bins  for  the  given  solar  flux  interval  (130  <  O  < 
130)  and  quiet  geomagnetic  conditions.  The  triangle  symbol  represents  bin  averages 
consisting  of  less  than  10  data  points. 
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resolution,  ensuring  we  used  only  quiet  time  measured  winds  (Kp  <  3),  and 
incorporating  two  additional  years  (1988-1989)  of  measured  FPI  data. 

5.2.2.3.  Solar  Cycle  Effects  (ISR  Measurements) 

In  addition  to  FPI  wind  data,  we  use  ISR-derived  neutral  winds  over  Arecibo. 
Although  these  data  are  not  the  primary  focus  of  this  study,  we  nevertheless  use  the  ISR 
results  to  supplement  the  Arecibo  FPI  data  analysis.  Although  it  is  possible  to  derive 
the  ISR  neutral  winds  between  altitudes  of  ~  200  -  400  km,  we  use  the  ISR  neutral  wind 
measured  at  the  calculated  airglow  emission  height,  ~  250  km,  to  compare  to  the  FPI 
measured  meridional  wind. 

The  ISR  neutral  winds  were  collected  during  January  1983  to  January  1987. 
Figure  49  shows  the  average  neutral  wind  response  to  two  levels  of  solar  activity  during 
the  equinox  and  June  solstice  periods.  This  data  set  contained  insufficient  data  in  the 
December  solstice  to  perform  a  valid  comparison  of  solar  activity  effects.  This  data  set 
is  characterized  by  low  to  moderate  solar  activity  conditions. 

Figure  49  shows  the  average  meridional  neutral  wind  for  low  (O  <  100)  and 
moderate  (3>  >  100)  solar  activity  during  the  equinox  and  June  solstice  periods  for  quiet 
geomagnetic  conditions.  Unlike  the  FPI  measurements,  the  ISR  is  not  restricted  to 
nighttime  measurements.  We  note  no  significant  variation  between  the  moderate  and 
low  flux  conditions  for  the  equinox  and  June  solstice  periods. 

5.2.2.4.  TIGCM93  Solar  Cycle  Activity  Model  Results 

Figures  50  and  5 1  show  the  Arecibo  TIGCM93  predicted  zonal  and  meridional 
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Figure  49.  ISR  neutral  winds  for  Arecibo,  Puerto  Rico,  for  various  seasons  for  low  and 
moderate  solar  activity.  The  average  neutral  wind  velocities  are  obtained  by  examining 
1/2  hour  time  bins  for  the  given  solar  flux  interval  (100  <  <I>  <  100)  and  quiet 
geomagnetic  conditions.  The  triangle  symbol  represents  bin  averages  consisting  of  less 
than  20  data  points. 
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neutral  winds  for  the  high,  moderate,  and  low  solar  activity  conditions.  The  model 
provides  results  for  the  December  solstice  (Day  =  360),  equinox  (Day  =  90),  and  June 
solstice  (Day  =  180)  periods.  Unlike  Arequipa,  the  combined  measured  (FPI)  and 
derived  (ISR)  wind  data  for  Arecibo  contain  a  sufficient  number  of  data  points  over  the 
solar  cycle  to  analyze  low  and  high  solar  flux  conditions.  We  use  the  relative  difference 
between  the  predicted  high  (O  =  195)  and  low  (^>  =  75)  solar  flux  conditions  to 
compare  to  the  experimental  wind  results. 

For  the  December  solstice.  Figure  50  shows  no  significant  variations  of  the 
TIGCM  zonal  wind  predictions  due  to  solar  activity  effects.  This  result  is  not  consistent 
with  the  experimental  results  since  the  measured  data  show  a  westward  perturbation 
from  low  to  high  solar  activity.  For  day  90,  the  model  predictions  show  no  wind 
variation  due  to  solar  activity,  a  result  which  agrees  with  the  experimental  data.  Last, 
the  model  results  for  day  180  show  a  westward  perturbation  from  high  to  low  solar 
activity  while  the  experimental  data  show  no  significant  variation  due  to  solar  cycle 
effects. 

Figure  51  shows  no  significant  variation  of  the  meridional  TIGCM93  wind 
predictions  to  solar  activity  except  for  day  180.  The  December  solstice  prediction 
shows  a  southward  wind  perturbation  from  solar  minimum  to  maximum  conditions. 
The  experimental  data  show  no  significant  wind  variation  due  to  solar  activity  for  all 
seasons.  Overall,  we  note  good  agreement  between  the  measured  and  TIGCM93 
predicted  response  of  the  Arecibo  winds  to  solar  cycle  effects  during  the  equinox  and 
June  solstice  periods  but  less  agreement  during  the  December  solstice  period. 
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Figure  50.  TIGCM93  eastward  wind  predictions  over  Arecibo,  Puerto  Rico,  for  various 
seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet  magnetic 
activity. 
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ARECIBO  TIGCM93 

0  0=195  A  0=  135  +0=  75  Kp  =  1 


Figure  51.  TIGCM93  northward  wind  predictions  over  Arecibo,  Puerto  Rico,  for 
various  seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet 
magnetic  activity. 
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5.2.2.5.  ErWM93  Solar  Cycle  Activity  Model  Results 

Figures  52  and  53  show  the  Arecibo  HWM93  modeled  zonal  and  meridional 
winds  for  the  high,  moderate,  and  low  solar  activity  conditions.  These  predicted  results 
are  for  the  December  solstice,  equinox,  and  June  solstice  periods.  We  use  the  predicted 
response  of  the  neutral  winds  during  high  (O  =  195)  and  low  (O  =  75)  solar  flux 
conditions  to  compare  to  the  experimental  results.  The  predicted  solar  cycle  effects  on 
the  zonal  winds  agree  with  our  experimental  data  except  for  the  postmidnight  period  for 
the  June  solstice  period.  The  HWM93  for  day  180  shows  a  westward  velocity 
perturbation  from  high  to  low  solar  activity,  a  result  not  consistent  with  the  measured 
wind  data  results.  A  possible  explanation  for  this  discrepancy  between  the  model  and 
observed  winds  is  the  HWM93  use  of  derived  ISR  winds  and  satellite  measured  winds 
in  addition  to  the  FPI  wind  measurements.  Therefore,  the  HWM93-predicted  westward 
perturbation  may  come  from  either  the  ISR  or  satellite  measurements,  or  a  combination 
of  the  ISR  and  satelHte  wind  measurements.  Figure  53  shows  the  meridional  predictions 
for  Arecibo.  These  data  show  no  significant  wind  variation  due  to  solar  activity.  This 
predicted  result  agrees  with  the  experimental  Arecibo  averaged  winds. 

5.3.  DATA  AND  MODEL  COMPARISONS 

Figures  54  -  57  show  comparisons  of  the  wind  results  from  the  TIGCM93,  the 
HWM93,  and  Arequipa  experimental  data.  We  compare  these  model  predictions  and 
experimental  average  winds  using  similar  solar  flux  and  magnetic  activity  conditions. 

Figure  54  compares  the  model  and  experimental  average  eastward  winds  over 
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ARECIBO  HWM93 

0  0=195  A  0=  135  +0=  75  Kp  =  2 


Figure  53.  HWM93  northward  wind  predictions  over  Arecibo,  Puerto  Rico,  for  various 
seasons  for  high,  moderate,  and  low  solar  activity  conditions  and  for  quiet  magnetic 
activity. 
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Arequipa  for  the  December  solstice,  equinox,  and  June  solstice  periods.  We  use  the 
model  winds  for  quiet  magnetic  activity  conditions,  Kp  =  1  for  the  TIGCM93  and  Kp  = 
2  for  the  HWM93,  and  moderate  solar  activity  except  during  the  June  solstice.  The 
Jiuie  solstice  experimental  wind  average  consists  of  wind  profiles  measured  during  high 
solar  activity;  hence,  we  use  the  modeled  winds  for  high  solar  flux  conditions  (Oth  = 
195).  As  expected,  the  HWM93  wind  predictions  agree  with  experimental  results  for  aU 
seasons.  The  TIGCM93  predictions  show  excellent  agreement  with  the  measured  winds 
in  the  June  solstice  and  equinox  periods.  However,  the  TIGCM93  shows  smaller 
eastward  winds  throughout  the  night  for  the  December  solstice  compared  to  the 
experimental  winds  during  the  same  time  period.  Figure  55  presents  the  modeled  and 
measured  meridional  winds  over  Arequipa.  The  model  predictions  agree  with  the 
experimental  wind  averages  during  the  equinox  and  December  solstice  periods,  but  the 
TIGCM93  predicts  slightly  larger  northward  winds  from  2300  LT  to  0230  LT 
compared  to  the  measured  data  during  the  June  solstice. 

Figures  56  and  57  show  the  model  and  experimental  winds  for  the  eastward  and 
northward  averages  over  Arecibo  for  the  December  solstice,  equinox,  and  June  solstice 
periods.  Where  possible,  we  used  model  results  with  similar  solar  flux  activity  to  that  of 
the  experimental  winds.  These  model  solar  flux  values  are  different  for  each  season  and 
are  provided  in  parentheses  below  the  average  experimental  solar  flux  and  magnetic 
activity  values. 

Figure  55  compares  the  TIGCM93  and  HWM93  eastward  predictions  to  the 
experimental  winds  over  Arecibo.  The  predicted  eastward  winds  agree  with  the 
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Figure  54.  Model  and  observed  eastward  wind  comparison  for  Arequipa,  Peru.  These 
comparisons  are  for  the  equinox  and  June  solstice.  The  TIGCM93  and  HWM93 
predictions  are  given  by  the  dashed  and  sohd  line,  respectively.  The  observed  wind 
results,  represented  by  the  open  circles  and  error  bars  (standard  deviation),  are  averaged 
FPI  data  for  moderate  solar  flux  conditions  and  quiet  magnetic  activity. 
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Figure  55.  Model  and  observed  northward  wind  comparison  for  Arequipa,  Peru.  These 
comparisons  are  for  the  equinox  and  June  solstice.  The  TIGCM93  and  HWM93 
predictions  are  given  by  the  dashed  and  solid  line,  respectively.  The  observed  wind 
results,  represented  by  the  open  circles  and  error  bars,  are  averaged  FPI  data  for 
moderate  solar  flux  conditions  and  quiet  magnetic  activity. 
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Figure  56.  Model  and  observed  eastward  \vind  comparison  for  various  seasons  for 
Arecibo,  Puerto  Rico.  The  TIGCM93  and  HWM93  predictions  are  given  by  the  dashed 
and  solid  line,  respectively.  The  Om  represents  the  solar  flux  value  used  for  the  model 
predictions.  The  average  wind  results,  represented  by  the  open  circles  and  error  bars, 
are  averaged  FPI  data  for  various  solar  flux  conditions  and  quiet  magnetic  activity. 
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Figure  57.  Model  and  observed  northward  wind  comparison  for  various  seasons  for 
Arecibo,  Puerto  Rico.  The  TIGCM93  and  HWM93  predictions  are  given  by  the  dashed 
and  solid  line,  respectively.  The  Oth  represents  the  solar  flux  value  used  for  the  model 
predictions.  The  average  wind  results,  represented  by  the  open  circles  and  error  bars, 
are  averaged  FPI  data  for  various  solar  flux  conditions  and  quiet  magnetic  activity. 
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measured  winds  during  the  premidnight  period  except  for  the  June  solstice  where  we 
note  poor  agreement  between  the  predicted  and  experimental  winds  throughout  the 
night.  We  also  note  that  the  TIGCM93  tends  to  predict  larger  eastward  winds  for  all 
seasons,  while  the  HWM93  shows  the  westward  reversal  occurring  at  an  earlier  time  for 
the  experimental  average  winds  during  the  equinox  and  June  solstice  periods.  Figure  57 
shows  good  agreement  between  the  Arecibo  model  northward  wind  predictions  and 
experimental  wind  averages. 
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CHAPTER  6 

MAGNETIC  ACTIVITY  EFFECTS  ON  THE  LOW  LATITUDE 
THERMOSPHERIC  WINDS 

6.1.  INTRODUCTION 

In  this  chapter  we  use  Fabry-Perot  measurements  at  Arequipa  (1988-1990)  and 
Arecibo  (1980-1989)  to  determine  the  magnetic  activity  effects  on  the  low  latitude  F- 
region  thermospheric  winds.  We  also  compare  the  magnetic  activity  response  of  the 
measured  winds  to  TIGCM93  and  HWM93  predictions. 

6.2.  MAGNETIC  ACnVITY  EFFECTS 


6.2.1.  Arequipa 

Figures  58  and  59  present  the  averaged  eastward  and  northward  neutral  winds 
over  Arequipa  during  magnetically  quiet  (Kp  <  3)  and  disturbed  (Kp  >  3)  conditions  for 
all  levels  of  solar  activity  during  the  equinox  and  June  solstice  periods.  The  Arequipa 
database  does  not  contain  enough  disturbed  wind  data  for  the  December  solstice  to 
adequately  determine  the  magnetic  activity  effects  on  the  neutral  winds  during  this 
period.  The  error  bars  represent  ±  1  standard  deviation  about  the  magnetically  quiet 
mean.  The  error  bars  for  the  disturbed  winds  are  not  shown  but  are  comparable  to 
those  of  the  quiet  wind.  The  results  in  Figure  58  show  a  westward  perturbation  when 
the  magnetic  activity  becomes  disturbed  for  the  equinox  and  June  solstice  periods.  The 
largest  decrease  in  eastward  winds  (~  25  -  50  m/s)  during  active  conditions  occurs  in  the 
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Figure  58.  Measured  eastward  wind  averages  over  Arequipa,  Peru,  for  the  equinox 
(March-April,  September-October)  and  June  solstice  (May-August)  periods  for  quiet 
(Kp  <  3)  and  disturbed  (Kp  >  3)  conditions.  The  solid  and  open  triangles  represent 
averages  based  on  less  than  10  data  points. 
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premidnight  period  for  both  equinox  and  June  solstice  seasons.  Figure  59  shows  the 
nighttime  average  neutral  meridional  winds  over  Arequipa  are  relatively  unchanged  for 
quiet  and  active  magnetic  periods  during  the  represented  seasons.  This  study  reports 
for  the  first  time  the  effects  of  geomagnetic  activity  on  the  ground-based  FPI  wind 
measurements  at  Arequipa,  Peru. 

Our  quiet  and  disturbed  zonal  wind  results  are  consistent  with  magnetic 
activity  studies  in  Wu  et  al.  [1994]  using  Dynamics  Explorer  2  (DE  2)  neutral  wind 
data.  The  DE  2  satellite  uses  data  collected  within  the  200  -  400  km  altitude  range 
during  1981  -  1983  for  day-  and  nighttime  conditions.  Although  most  DE  2  data 
were  collected  during  the  December  solstice,  some  DE  2  data  were  collected  during 
September  -  October,  This  data  set  allows  us  to  compare  DE  2  results  with  our  FPI 
equinox  winds.  The  DE  2  data  were  subdivided  into  geomagnetically  active  (BCp  > 
3)  and  quiet  (Kp  <  3)  periods  to  analyze  the  response  of  the  neutral  wind  to 
magnetic  activity  effects  [Wu  et  al,  1994].  Wu  et  al.  [1994]  showed  a  westward 
perturbation  when  the  magnetic  activity  became  disturbed,  a  result  that  agrees  with 
our  FPI  disturbed  observations.  The  DE  2  data  contained  too  few  nighttime 
disturbed  meridional  wind  data  points  to  properly  investigate  magnetic  activity 
effects  on  the  northward  winds. 

6.2.2.  Arecibo 

We  have  also  averaged  the  eastward  and  northward  neutral  winds  for  quiet  and 
disturbed  conditions  over  Arecibo  for  each  season.  Figure  60  shows  the  averaged 
eastward  winds  for  all  levels  of  solar  activity  during  the  December  solstice,  equinox,  and 
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Figure  59.  Measured  northward  wind  averages  over  Arequipa,  Peru,  for  the  equinox 
(March-April,  September-October)  and  June  solstice  (May-August)  periods  for  quiet 
(I^  <  3)  and  disturbed  (Kp  >  3)  conditions.  The  solid  and  open  triangles  represent 
averages  based  on  less  than  10  data  points. 
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Figure  60.  Measured  eastward  wind  averages  over  Aiecibo,  Puerto  Rico,  for  various 
seasons  for  quiet  (Kp  <  3)  and  disturbed  (Kp  >  3)  conditions.  The  solid  and  open 
triangles  represent  averages  based  on  less  than  10  data  points. 
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June  solstice  periods.  As  the  magnetic  activity  becomes  disturbed,  the  data  show  a 
westward  perturbation  during  all  nighttime  hours  for  the  December  solstice  and  equinox 
periods  while  we  note  a  westward  permrbations  in  the  postmidnight  period  during  the 
June  solstice.  Figure  61  shows  no  significant  variation  between  the  quiet  and  disturbed 
northward  winds. 

Burnside  et  al.  [1991]  showed  westward  perturbations  over  Arecibo  during 
storm-time  conditions.  This  previous  study  agrees  with  our  measured  wind  results. 
However,  these  studies  concentrated  on  a  few  disturbed  days  and  did  not  provide  any 
indication  of  the  average  neutral  wind  response  to  magnetic  activity  effects  at  Arecibo. 

Buonsanto  and  Foster  [1993]  and  Burnside  et  al.  [1991]  analyzed  the 
meridional  wind  response  to  disturbed  conditions  over  Arecibo.  Although  these  studies 
showed  that  the  low  latitude  meridional  wind  responds  to  storm  time  effects,  they 
showed  both  northward  and  southward  perturbations  disturbed  conditions.  Even 
though  these  studies  concentrated  on  a  few  disturbed  days,  they  are  consistent  with  our 
results.  These  northward  and  southward  wind  perturbations  would  tend  to  cancel  each 
other  over  a  long-term  average.  This  is  consistent  with  our  results  over  Arecibo 
because  we  found  no  significant  variation  in  the  average  meridional  wind  response  to 
magnetic  activity  effects. 

In  addition  to  FPI  wind  data,  we  use  ISR-derived  neutral  winds  over  Arecibo. 
The  ISR  neutral  winds  were  collected  during  January  1983  to  January  1987.  Figure  62 
shows  the  quiet  and  disturbed  wind  averages  for  the  December  solstice,  equinox,  and 
June  solstice.  We  note  no  significant  variation  between  the  December  solstice  disturbed 
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Figure  61.  Measured  northward  wind  averages  over  Arecibo,  Puerto  Rico,  for  various 
seasons  for  quiet  (Kp  <  3)  and  disturbed  (Kp  >  3)  conditions.  The  solid  and  open 
triangles  represent  averages  based  on  less  than  10  data  points. 
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Figure  62.  Measured  ISR  neutral  winds  over  Arecibo,  Puerto  Rico,  for  various  seasons 
for  quiet  (Kp  <  3)  and  disturbed  (Kp  >  3)  conditions.  The  solid  and  open  triangles 
represent  averages  based  on  less  than  40  data  points. 
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and  quiet  winds  from  1600  LT  to  0900  LT,  but  the  data  show  a  southward  perturbation 
from  quiet  to  disturbed  magnetic  activity  from  1000  LT  to  1500  LT.  The  data  also 
show  no  significant  magnetic  activity  effects  during  the  equinox  and  June  periods.  The 
ISR  and  FPI  results  indicate  the  response  of  the  average  meridional  winds  is  not 
affected  by  magnetic  activity. 

6.3.  SHORT-TERM  AND  EXTENDED  MAGNETIC  ACTIVITY 

We  analyze  the  short-term  (3  hour)  and  extended  time  (12  hour)  magnetic 
activity  to  fully  investigate  the  neutral  wind  response  to  storm-time  events.  We  define 
three  magnetic  activity  cases,  which  cover  all  storm-time  scenarios.  These  magnetic 
activity  classifications  are  listed  in  TABLE  16. 

The  first  magnetic  activity  classification  is  the  extended  quiet  condition.  This 


TABLE  16.  Magnetic  Activity  Classification 


Classification  Name 

Kp  Index 

Magnetic  Activity  Name 

mQ 

12hrKp<3 

Kp<3 

extended  quiet 

2:q,d 

12hrKp<3 

Kp>3 

short-term  disturbed 

5:d,d 

12hrKp>3 

Kp>3 

extended  disturbed 
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case  is  defined  by  a  magnetically  quiet  12-hour  Kp  and  3-hour  Kp.  The  12-hour  Kp  is  the 
average  of  the  three  consecutive  Kp  values  preceding  the  Kp  value  at  the  time  of  the 
wind  measurement.  The  3-hour  Kp  is  the  magnetic  activity  defined  at  the  time  of  wind 
measurement.  The  second  classification  highlights  short-term  disturbed  effects.  This 
case  is  defined  by  the  previous  magnetic  activity,  which  is  quiet  (12-hour  Kp  <  3)  and 
then  becomes  disturbed  (Kp  >  3).  The  last  case  defines  the  condition  where  the  previous 
and  short-term  magnetic  activity  are  disturbed.  This  is  the  first  time  this  type  of 
magnetic  activity  classification  has  been  applied  to  investigate  the  response  of  the 
thermospheric  neutral  winds  to  magnetic  activity  effects. 

The  extended  quiet  and  disturbed  classifications  define  the  extreme  range  of  the 
magnetic  activity  conditions.  If  no  wind  variation  is  noted  between  the  extended  quiet 
and  disturbed  cases,  we  conclude  there  will  be  no  short-term  disturbed  effects  since  this 
case  is  included  in  the  disturbed  baseline. 

6.4.  NEUTRAL  WIND  RESPONSE  TO  MAGNETIC  ACTIVITY 

6.4.1.  Arequipa  and  Arecibo  Northward  Wind 

Figures  63  and  64  present  the  extended  quiet  and  disturbed  northward  winds 
over  Arequipa  and  Arecibo  for  all  levels  of  solar  activity.  Figure  63  shows  the  results 
for  the  June  solstice  and  Figure  64  shows  the  results  for  the  equinox  period.  Our  data 
set  does  not  contain  enough  statistically  reliable  averages  during  the  December  solstice 
to  conduct  a  proper  analysis  of  magnetic  activity  effects  for  both  Arequipa  and  Arecibo. 
These  graphs  show  the  average  wind  profile  and  spline  fits  for  the  average  solar  flux 
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Figure  63.  June  solstice  extended  quiet  (£Q,Q)  and  extended  disturbed  (2D,D) 
northward  wind  comparisons  for  Aiequipa  and  Aiecibo.  The  solid  line  is  the  spline  fit 
for  the  ZQ,Q  case  and  the  dashed  line  is  the  spline  fit  for  the  ED,D  case.  The  open  and 
solid  triangles  represent  averages  consisting  of  less  than  10  wind  data  points. 
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Figure  64.  Equinox  extended  quiet  (£Q,Q)  and  extended  disturbed  (XD,D)  northward 
wind  comparisons  for  Arequipa  and  Arecibo.  The  solid  line  is  the  spline  fit  for  the 
XQ,Q  case  and  the  dashed  line  is  the  spline  fit  for  the  XD,D  case.  The  open  and  solid 
triangles  represent  averages  consisting  of  less  than  10  wind  data  points. 
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(<I>),  average  short-term  magnetic  activity  (Kp),  and  average  extended  magnetic  activity 
(12hrKp)  conditions  for  Arequipa  and  Arecibo.  The  open  symbols  represent  the 

extended  quiet  average  wind  profile,  and  the  filled  symbols  represent  the  extended 
disturbed  average  wind  profile.  The  solid  line  is  the  spline  fit  for  the  quiet  average 
profile,  and  the  dashed  line  is  the  spline  fit  for  the  disturbed  average  profile.  The  open 
and  filled  triangle  symbols  represent  average  wind  values  consisting  of  less  than  ten  data 
points. 

Figures  63  and  64  show  no  significant  variation  between  the  extended  quiet  and 
disturbed  northward  averages.  Because  these  comparisons  were  made  using  the 
extreme  magnetic  activity  cases,  we  conclude  that  the  average  northward  wind  response 
for  Arequipa  and  Arecibo  is  not  affected  by  magnetic  activity. 

6.4.2.  Arequipa  and  Arecibo  Eastward  Wind 

Figures  65  and  66  present  the  extended  quiet  and  disturbed  eastward  winds 
over  Arequipa  and  Arecibo.  We  use  the  same  graph  setups  described  in  the  previous 
section  but  applied  to  the  eastward  quiet  and  disturbed  wind  averages.  The  solid  and 
the  dashed  lines  are  the  spline  fits  for  the  quiet  and  disturbed  average  winds, 
respectively. 

The  near- June  solstice  and  equinox  periods  show  a  significant  westward  velocity 
perturbation  for  the  disturbed  condition  for  Arequipa  and  Arecibo.  We  compared  the 
extended  quiet  to  the  short-term  disturbed  condition  for  the  June  solstice  periods  in 
Figure  67  to  further  isolate  the  magnetic  activity  effects  on  the  averaged  eastward 
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May-August  Zonal  Wind  Averages  All  Fluxes 

Extended  Quiet  Extended  Disturbed 

O  -  2Q,Q  vs  m-  1D,D 


Local  Time 


Figure  65.  June  solstice  extended  quiet  (2Q,Q)  and  extended  disturbed  (ZD,D) 
eastward  wind  comparisons  for  Arequipa  and  Arecibo.  The  solid  line  is  the  spline  fit  for 
the  2Q,Q  case  and  the  dashed  line  is  the  spline  fit  for  the  ED,D  case.  The  open  and 
solid  triangles  represent  averages  consisting  of  less  than  10  wind  data  points. 
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Equinox  Zonal  Wind  Averages  All  Fluxes 

Extended  Quiet  Extended  Disturbed 

O  -  IQ,Q  vs  •  -  ZD,D 


Local  Time 


Figure  66.  Equinox  extended  quiet  (EQ,Q)  and  extended  disturbed  (LD,D)  eastward 
wind  comparisons  for  Arequipa  and  Arecibo.  The  solid  line  is  the  spline  fit  for  the 
ZQ,Q  case  and  the  dashed  line  is  the  spline  fit  for  the  XD,D  case.  The  open  and  solid 
triangles  represent  averages  consisting  of  less  than  10  wind  data  points. 
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winds.  The  wind  databases  for  Arequipa  and  Arecibo  do  not  contain  enough  short-term 
disturbed  wind  data  to  perform  a  valid  comparison  of  magnetic  activity  effects  on  the 
winds  for  the  equinox  period. 

Figure  67  shows  no  significant  variation  between  the  extended  quiet  and  short¬ 
term  disturbed  eastward  winds  over  Arequipa.  For  Arecibo,  we  note  a  small  westward 
velocity  perturbation  as  the  short-term  becomes  disturbed  in  the  postmidnight  period 
during  the  June  solstice. 

6.5  MAGNETIC  ACTIVITY  CUBIC  SPLINE  AND 
MODEL  PREDICTIONS 


6.5.1.  Arequipa 

6.5. 1. 1 .  Cubic  Spine  Predictions 

We  now  investigate  these  same  magnetic  activity  conditions  using  the  cubic 
spline  fits.  Figures  68  and  69  show  the  spline  fits  for  the  northward  and  eastward 
winds  over  Arequipa  for  the  extended  quiet  (ZQ,Q),  short-term  disturbed  (EQ,D),  and 
extended  disturbed  (ZD,D)  magnetic  activity  cases.  The  solar  flux  is  set  to  moderate 
flux  conditions,  O  =  150,  to  eliminate  variations  due  to  solar  activity. 

Figure  68  shows  the  cubic  spline  predictions  for  the  northward  winds  during  the 
equinox  (day  90)  and  December  solstice  (day  180)  periods.  The  June  solstice  (day  360) 
prediction  is  not  shown  due  to  the  limited  number  of  wind  data  points  for  this  season. 
These  spline  predictions  show  no  significant  variations  in  the  neutral  wind  response  due 
to  magnetic  activity  for  days  90  and  180. 
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AREQUIPA  Spline  Fit 


OIQ,Q  12hrKp  =  2,  Kp  =  2 
Am,D  12hrKp  =  2,  Kp  =  5 
+  ED,D  12hrKp  =  5,  Kp  =  5  0=  150 


Figure  67.  June  solstice  extended  quiet  (ZQ,Q)  and  short-term  disturbed  (EQ,D) 
eastward  wind  comparisons  for  Arequipa  and  Arecibo.  The  solid  line  is  the  spline  fit  for 
the  XQ,Q  case  and  the  dashed  line  is  the  spline  fit  for  the  XQ,D  case.  The  open  and  solid 
triangles  represent  averages  consisting  of  less  than  10  wind  data  points. 


Northward  Winds  (m/s) 
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ARECIBO  Spline  Fit 

O  ZQ.Q  12hrKp  =  2,  Kp  =  2 
A1Q,D  12hrKp  =  2,  Kp  =  5 
+  ED,D  12hrKp  =  5,  Kp  =  5  0=150 


Figure  68.  Cubic  spline  predictions  for  the  northward  winds  over  Arequipa  for  equinox 
(day  =  90)  and  June  solstice  (day  =  180)  for  moderate  solar  flux  (^  =  150)  conditions. 
The  spline  fit  results  are  for  the  extended  quiet  (LQ,Q),  short-term  disturbed  (XQ,D), 
and  the  extended  disturbed  case  (ZD,D). 
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Figure  69  shows  the  cubic  spline  predictions  for  the  eastward  wind  over 
Arequipa.  We  calculated  spline  fits  for  days  90  and  180.  Day  360  results  are  not  shown 
due  to  the  limited  number  of  wind  data  points  collected  during  this  period.  For  days  90 
and  180,  the  data  show  no  significant  variation  between  the  extended  quiet  and  the 
short-term  disturbed  cases;  however,  the  spline  prediction  shows  a  westward 
perturbation  through  the  night  when  the  extended  magnetic  activity  becomes  disturbed. 

6.5. 1.2.  TIGCM93  Predictions 

Figures  70  and  71  present  the  seasonal  Arequipa  TIGCM93  eastward  and 
northward  predictions  for  magnetically  quiet  and  disturbed  conditions  defined  by  Kp  =  1 
and  Kp  =  6,  respectively.  We  use  the  predicted  results  for  high  solar  activity  (^>  =  195). 
For  each  TIGCM93  simulation,  the  model  is  run  to  diumally  reproducible  solutions  or 
steady-state  conditions.  The  steady-state  condition  means  that  the  value  of  a  modeled 
field  (e.g.,  temperature,  neutral  wind,  electron  density,  etc.)  at  a  given  TIGCM93  grid 
point  at  a  given  universal  time  falls  within  a  few  percent  of  the  value  calculated  on  the 
previous  day.  Since  these  model  runs  were  produced  using  steady-state  conditions,  the 
TIGCM93-modeled  quiet  and  disturbed  winds  are  most  like  the  extended  quiet  and 
disturbed  magnetic  activity  outline  in  TABLE  16.  In  chapter  4  we  found  the  average 
airglow  emission  height  to  be  near  250  km  through  the  night;  therefore,  we  use  the 
predicted  winds  at  250  km. 

Similar  to  the  approach  used  for  the  experimental  wind,  we  analyze  the 
predicted  wind  response  to  magnetic  activity  effects  during  the  equinox  and  June 
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AREQUIPA  Spline  Fit 

O  IQ.Q  12hrKp  =  2,  Kp  =  2 
A  EQ,D  12hrKp  =  2,  Kp  =  5 
+  1D,D  12hrKp  =  5,  Kp  =  5  150 


Figure  69.  Cubic  spline  predictions  for  the  eastward  winds  over  Arequipa  for  equinox 
(day  =  90)  and  June  solstice  (day  =  180)  for  moderate  solar  flux  (O  =  150)  conditions. 
The  spline  fit  results  are  for  the  extended  quiet  (XQ.Q),  short-term  disturbed  (XQ,D), 
and  the  extended  disturbed  case  (ZD,D). 
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solstice  period.  We  focus  on  the  nighttime  TIGCM93  wind  predictions  because  we 
want  to  compare  the  model  results  to  the  FPI  wind  measurements. 

Figure  70  presents  the  eastward  wind  predictions.  We  note  a  westward  velocity 
perturbation  for  most  of  the  night  for  the  disturbed  wind  during  aU  seasons.  We  see  no 
significant  variation  between  the  quiet  and  disturbed  wind  from  0300  LT  to  the  early 
morning  hours.  For  the  northward  wind  prediction  in  Figure  71,  we  note  no  significant 
variation  between  the  quiet  and  disturbed  wind  for  days  90  and  180,  but  we  do  note  a 
southward  velocity  perturbation  as  the  magnetic  activity  becomes  disturbed  in  the 
postmidnight  period  for  day  360. 

6.5. 1.3  HWM93  Predictions 

Figures  72  and  73  present  the  seasonal  HWM93  model  predictions  for 
magnetically  quiet  and  disturbed  conditions  defined  by  Kp  =  2  and  Kp  =  5,  respectively. 
We  use  predicted  results  for  moderate  solar  activity  (^>  =  150).  The  HWM93  uses  the 
daily  Ap  magnetic  activity  index;  therefore,  this  model  emphasizes  the  long-term 
magnetic  activity  effects. 

The  eastward  predictions  in  Figure  72  show  no  significant  variations  between 
the  quiet  and  disturbed  winds  in  the  premidnight  period  for  all  seasons;  however,  the 
model  winds  show  a  small  westward  velocity  perturbation  as  the  magnetic  activity 
becomes  disturbed  in  the  postmidnight  period  for  all  seasons.  The  Arequipa  northward 
predictions  in  Figure  73  show  no  significant  variation  between  the  quiet  and  disturbed 
wind  throughout  the  night  for  all  seasons. 
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AREQUIPA  TIGCM93 

O  Kp  =  1  A  Kp  =  6  0=  195 


Figure  70.  TIGCM93  quiet  (Kp  =  1)  and  disturbed  (Kp  =  6)  eastward  wind  predictions 
over  Arequipa  for  all  seasons.  These  predictions  are  for  high  solar  flux  (<I>  =  195) 
conditions. 


Northward  Winds  (m/s) 
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AREQUIPA  TIGCM93 

O  Kp  =  1  A  Kp  =  6  0=  195 


Figure  71.  TIGCM93  quiet  (Kp  =1)  and  disturbed  (Kp  =  6)  northward  wind  predictions 
over  Arequipa  for  all  seasons.  These  predictions  are  for  high  solar  flux  (O  =  195) 
conditions. 


Eastward  Winds  (m/s) 
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AREQUIPA  HWM93 

O  Kp  =  2  A  Kp  =  5  0=  150 


Figure  72.  HWM93  quiet  (Kp  =  2)  and  disturbed  (Kp  =  5)  eastward  wind  predictions 
over  Arequipa  for  all  seasons.  These  predictions  are  for  moderate  solar  flux  (O  =  150) 
conditions. 


155 


AREQUIPA  HWM93 

O  Kp  =  2  A  Kp  =  5  0=  150 


Figure  73.  HWM93  quiet  (Kp  =  2)  and  disturbed  (Kp  =  5)  northward  wind  predictions 
over  Arequipa  for  all  seasons.  These  predictions  are  for  moderate  solar  flux  (C>  =  150) 
conditions. 
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6.5.2.  Arecibo 

6.5.2. 1.  Cubic  Spine  Predictions 

Figure  74  shows  the  spline  fits  for  northward  wind  over  Arecibo  during  the 
December  solstice  (day  360),  equinox  (day  90),  and  June  solstice  (day  180).  We  do  not 
show  the  short-term  disturbed  (21Q,D)  effects  for  days  360  and  180  due  to  the  limited 
short-term  disturbed  wind  data  for  the  December  and  June  solstices.  These  spline  fits 
do  not  show  noticeable  magnetic  activity  effects  on  the  northward  component  for  any 
season.  The  cubic  spline  fits  for  the  northward  wind  agree  with  the  experimental  results 
under  similar  magnetic  activity  conditions. 

Figure  75  shows  the  cubic  spline  fits  for  the  zonal  winds  for  the  December 
solstice,  equinox,  and  June  solstice.  In  contrast  to  Arequipa,  the  Arecibo  spline  results 
in  Figure  75  show  a  westward  velocity  perturbation  when  the  short-term  magnetic 
activity  becomes  disturbed  for  days  90  and  180.  This  westward  perturbation  continues 
with  extended  disturbed  magnetic  activity  conditions.  For  day  360,  we  note  no 
significant  variation  in  the  spline  result  due  to  magnetic  activity  effects. 

This  westward  perturbation  is  consistent  with  theoretical  results  by  Blanc  and 
Richmond  [1980].  Using  a  first  principal  theoretical  model,  they  showed  the 
thermospheric  winds  above  120  km  produced  by  auroral  heating  during  magnetic  storms 
are  driven  equatorward.  The  transport  of  angular  momentum  by  these  winds  produces  a 
westward  motion  with  respect  to  the  earth.  Hence,  magnetic  activity  effects  reach 


Arecibo  earlier  than  Arequipa. 


Northward  Winds  (m/s) 
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ARECIBO  Spline  Fit 

O  m,Q  12hrKp  =  2,  Kp  =  2 
ATQ,D  12hrKp  =  2,  Kp  =  5 
+  XD,D  12hrKp  =  5,  Kp  =  5  0=  150 


Figure  74.  Cubic  spline  predictions  for  the  northward  winds  over  Arecibo  for  all 
seasons  and  moderate  solar  flux  (<>  =  150)  conditions.  The  spline  fit  results  are  for  the 
extended  quiet  (2Q,Q),  short-term  disturbed  (2Q,D),  and  the  extended  disturbed 
(XD,D)  cases. 


Eastward  Winds  (m/s) 
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ARECIBO  Spline  Fit 

O  1Q,Q  12hrKp  =  2.  Kp  =  2 
A  EQ,D  12hrKp  =  2,  Kp  =  5 
+  ED.D  12hrKp  =  5.  Kp  =  5  0=  150 


19  21  23  01  03  05  07 

Local  Time 


Figure  75.  Cubic  spline  predictions  for  the  eastward  winds  over  Arecibo  for  all  seasons 
and  moderate  solar  flux  (O  =150)  conditions.  The  spline  fit  results  are  for  the  extended 
quiet  (EQ,Q),  short-term  disturbed  (IQ,D),  and  the  extended  disturbed  case  (SD,D). 
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6.5.22.  TIGCM93  Predictions 

Figures  76  and  77  show  the  Arecibo  eastward  and  northward  TIGCM93 
predicted  winds  for  quiet  and  disturbed  magnetic  activity  and  low  solar  activity.  The 
Arecibo  eastward  winds  in  Figure  76  show  a  significant  westward  velocity  perturbation 
for  the  disturbed  wind  for  most  of  the  night.  We  note  no  significant  variation  between 
the  quiet  and  disturbed  winds  from  0400  LT  to  early  morning.  The  northward  wind 
predictions  in  Figure  77  show  a  southward  velocity  perturbation  through  the  night  for 
days  360  and  90,  while  we  note  a  southward  velocity  perturbation  until  0500  LT  for  day 
180. 

6.5.2.3.  HWM93  Predictions 

Figures  78  and  79  show  the  HWM93  predictions  for  the  eastward  and 
northward  winds  over  Arecibo.  Figure  78  presents  the  predicted  eastward  results  for 
quiet  and  disturbed  winds  for  moderate  solar  activity  conditions.  We  note  a  small 
westward  velocity  perturbation  for  most  of  the  night  during  all  seasons.  The  relative 
velocity  difference  between  the  quiet  and  disturbed  eastward  winds  is  ~  25  m/s  for  all 
seasons.  The  Arecibo  northward  predictions  in  Figure  79  show  no  significant  variation 
between  the  quiet  and  disturbed  wind  throughout  the  night  for  any  season. 
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ARECIBO  TIGCM93 

O  Kp  =  1  A  Kp  =  6  0=  75 


Figure  76.  TIGCM93  quiet  (Kp  =  1)  and  disturbed  (Kp  =  6)  eastward  wind  predictions 
over  Arecibo  for  all  seasons.  These  predictions  are  for  low  solar  flux  (O  =  75) 
conditions. 


Northward  Winds  (m/s) 


161 


ARECIBO  TIGCM93 

O  Kp  =  1  A  Kp  =  6  0=  75 


Figure  77.  TIGCM93  quiet  (Kp  =  1)  and  disturbed  (Kp  =  6)  northward  wind  predictions 
over  Arecibo  for  all  seasons.  These  predictions  are  for  low  solar  flux  (O  =  75) 
conditions. 
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ARECIBO  HWM93 

O  Kp  =  2  A  Kp  =  5  0=  150 


Day  =  360 


Local  Time 


Figure  78.  HWM93  quiet  (Kp  =  2)  and  disturbed  (Kp  =  5)  eastward  wind  predictions 
over  Arecibo  for  all  seasons.  These  predictions  are  for  moderate  solar  flux  (<b  =  150) 
conditions. 
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ARECIBO  HWM93 

O  Kp  =  2  A  Kp  =  5  0=  150 


Figure  79.  HWM93  quiet  (Kp  =  2)  and  disturbed  (Kp  =  5)  northward  wind  predictions 
over  Arecibo  for  all  seasons.  These  predictions  are  for  moderate  solar  flux  (^>  =  150) 
conditions. 
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CHAPTER? 

CONCLUSIONS 

7.1.  CONCLUSIONS 

The  primary  purpose  of  this  thesis  is  to  investigate  magnetic  activity  effects  on 
the  low  latitude  thermospheric  neutral  winds.  To  properly  investigate  storm-time 
effects  on  the  neutral  winds,  we  first  calculated  the  average  airglow  emission  height  for 
different  solar  cycle  and  magnetic  activity  conditions.  Next,  we  analyzed  the  solar  cycle 
effects  on  the  seasonal  averages.  The  following  paragraphs  outline  our  conclusions  on 
the  airglow  emission  heights,  solar  cycle  effects,  and  magnetic  activity  analysis. 

7.1.1.  Airglow  Emission  Heights 

We  used  electron  density  measurements  at  Arecibo,  Puerto  Rico,  and  Jicamarca, 
Peru,  MSISE90  model  predictions,  and  algorithms  developed  by  Link  and  Cogger 
[1988]  to  calculate  the  average  airglow  emission  height  over  the  low  latitude 
ionosphere.  This  calculated  emission  height  provides  us  the  average  altitude  of  the 
thermospheric  wind  measurement  and  the  reference  height  for  modeled  and  measured 
wind  comparisons.  Our  study  shows  the  average  airglow  emission  height  does  not  vary 
with  magnetic  activity.  The  airglow  height  above  Arequipa  does  not  change  with  solar 
activity;  however,  we  note  lower  average  heights  from  solar  maximum  to  minimum  over 
Arecibo.  Overall,  the  average  airglow  emission  height  is  ~  250  km.  Our  airglow 
emission  height  calculations  are  consistent  with  observations  by  Barrage  et  al.  [1990]  in 
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which  he  calculated  the  December  solstice  airglow  emission  peak  to  be  ~  270  km  at 
0100  LT. 

7.1.2.  Solar  Activity  Effects 

The  averaged  meridional  wind  is  not  affected  by  solar  cycle  effects  at  either 
Arequipa  or  Arecibo.  The  June  solstice  at  Arequipa  shows  an  eastward  wind  increase 
of  ~  50  m/s  during  the  postmidnight  period  from  moderate  to  high  solar  flux  conditions. 
The  Arecibo  December  solstice  period  shows  an  eastward  decrease  of  ~  50  m/s  during 
the  premidnight  period  from  solar  minimum  to  maximum  conditions. 

Predicted  TIGCM93  meridional  wind  response  to  solar  activity  agrees  with  the 
experimental  results  at  Arequipa  and  Arecibo.  However,  the  TIGCM93  predicts  higher 
eastward  winds  in  the  premidnight  period  from  low  to  high  solar  flux  conditions  at 
Arequipa,  and  the  modeled  eastward  wind  response  to  solar  cycle  effects  over  Arecibo 
does  not  agree  with  our  experimental  results. 

Predicted  HWM93  meridional  wind  response  to  solar  activity  agrees  with  the 
experimental  results  at  Arequipa  and  Arecibo.  However,  the  model  predicts  increased 
eastward  winds  (~  25  m/s)  in  the  postmidnight  period  from  low  to  high  solar  flux 
condition  at  Arequipa  during  the  equinox  and  June  solstice  periods.  These  predictions 
do  not  agree  with  the  experimental  results  for  the  equinox  period.  The  model  results  for 
the  December  solstice  show  a  smaller  wind  response  (  ~  25  m/s  )  between  the  moderate 
to  high  solar  activity  average  winds  than  that  shown  in  our  experimental  analysis  (~  50 
m/s).  Over  Arecibo,  the  zonal  wind  predictions  agree  with  the  December  solstice 
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measured  winds  but  not  with  the  June  solstice  period.  The  Arecibo  HWM93  June 
predictions  show  an  eastward  increase  of  ~  50  m/s  in  the  postmidnight  period  from  low 
to  high  solar  flux  conditions  while  our  June  solstice  experimental  FPI  data  show  no 
variation  due  to  solar  cycle  effects. 

Our  results  of  the  response  of  the  meridional  neutral  winds  over  Arequipa  and 
Arecibo  to  solar  cycle  effects  are  consistent  with  previous  wind  studies  conducted  at 
these  locations  [Biondi  et  al.,  1991;  Burnside  and  Tepley,  1989].  Our  FPI-measured 
data  showed  solar  cycle  effects  on  the  zonal  winds  at  Arequipa  and  Arecibo.  Previous 
studies  outlined  in  Biondi  et  al.  [1991]  showed  larger  eastward  winds  from  solar 
minimum  to  maximum  conditions,  a  result  consistent  with  our  solar  cycle  analysis.  For 
Arecibo,  Burnside  and  Tepley  [1989]  showed  the  average  zonal  wind  is  not  affected  by 
solar  cycle  effects,  a  result  that  does  not  agree  with  our  study.  We  enhanced  the  solar 
cycle  study  performed  in  Burnside  and  Tepley  [1989]  by  using  only  quiet- time  (Kp  <  3) 
wind  data,  incorporating  two  additional  years  (1988  and  1989)  of  FPI  wind  data,  and 
implementing  a  finer  time  resolution  (1/2-hour  time  bins). 

7.1.3.  Magnetic  Activity  Effects 

The  averaged  meridional  wind  is  not  affected  by  magnetic  activity  at  either 
Arequipa  or  Arecibo  based  on  the  Kp  index.  We  note  that  a  westward  perturbation  or 
eastward  decrease  of  the  average  zonal  wind  due  to  magnetic  activity  occurs  on  a  short 
time  scale  at  Arecibo  but  not  at  Arequipa.  We  conclude  that  magnetic  activity  effects 
reach  Arecibo  before  Arequipa.  Both  locations  show  large  westward  perturbations  for 
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long-term  disturbed  effects.  This  westward  wind  response  to  disturbed  magnetic 
activity  conditions  has  been  shown  in  previous  experimental  studies  [Burnside  et  al, 
1991;  Wu  etal.,  1994]  and  theoretical  simulations  [Blanc  and  Richmond,  1980], 

The  predicted  TIGCM93  zonal  wind  response  to  magnetic  activity  agrees  with 
experimental  results  in  the  premidnight  period  at  both  locations.  The  modeled 
meridional  response  to  magnetic  activity  agrees  with  the  experimental  results  at 
Arequipa  but  not  at  Arecibo.  The  Arecibo  TIGCM93  meridional  predictions  show  a 
northward  decrease  of  ~  40  -  50  m/s  between  the  quiet  and  disturbed  winds  in  the 
premidnight  period  for  all  seasons. 

The  predicted  HWM93  zonal  wind  response  to  magnetic  activity  agrees  with 
experimental  results  in  the  postmidnight  period  for  Arequipa  and  all  nighttime  hours  at 
Arecibo;  however,  the  model  underestimates  the  disturbed  wind  response.  The 
predicted  meridional  wind  response  agrees  with  the  experimental  results. 

7.2.  FUTURE  WORK 

Clearly,  more  neutral  wind  data  from  FPIs,  ISRs,  and  satellites  are  needed  to 
fuUy  investigate  the  solar  cycle  and  magnetic  activity  effects.  A  few  ideas  for  future 
research  are  presented  below. 

This  thesis  presented  the  climatological  averages  of  the  neutral  wind  response  to 
disturbed  magnetic  activity  conditions.  We  compared  these  averages  to  the  HWM93 
and  'nGCM93  model  predictions  to  find  out  how  well  these  models  could  reproduce 
the  experimental  results.  One  possible  research  opportunity  would  be  to  vary  the  inputs 


168 


(tides,  solar  flux,  geomagnetic  activity,  etc.)  to  the  first  principal  model,  the  TIGCM93, 
to  fully  investigate  the  thermospheric  physics  and  coupling  mechanisms.  This 
investigation  may  provide  further  clues  as  to  why  the  average  experimental  meridional 
wind  response  does  not  vary  between  quiet  and  disturbed  conditions.  A  more  complete 
analysis  would  require  other  first  principal  models  currently  used  in  the  space  science 
community. 

Recent  investigations  by  Fejer  [1993]  described  the  coupling  between  the  F- 
region  nighttime  zonal  neutral  winds  and  plasma  drifts  over  Arecibo  during  magnetically 
quiet  conditions.  Besides  the  experimental  results,  Blanc  and  Richmond  [1980] 
describe  model  predictions  of  the  neutral  winds  and  plasma  drifts  during  magnetically 
disturbed  conditions.  A  recent  paper  by  Coley  et  al.  [1994]  uses  measurements  of  the 
zonal  neutral  wind  and  zonal  ion  drifts  by  the  DE  2  satellite  to  conclude  that  on  average 
the  zonal  ion  drift  speed  exceeds  the  zonal  neutral  wind  speed.  These  DE  2  data  were 
collected  during  April  1981  through  February  1983,  a  period  characterized  by  high  solar 
flux  conditions.  Tlie  data  were  collected  in  the  300  -  400  km  altitude  range  near  the  dip 
equator.  During  2000  -  0500  LT,  Coley  et  al.  [1994]  showed  the  ion  drift  is  more 
strongly  eastward  than  the  neutral  wind.  This  result  is  not  consistent  with  the 
experimental  results  in  Fejer  [1993]  and  model  predictions  by  Blanc  and  Richmond 
[1980].  This  discrepancy  could  be  further  explored  using  the  FPI  measurements  at 
Arecibo  and  Arequipa  and  zonal  drifts  measured  by  the  ISRs  at  Arecibo  and  Jicamarca. 
A  detailed  comparison  of  the  neutral  winds  and  zonal  drifts  for  quiet  and  disturbed 
conditions  would  be  useful  for  investigating  the  results  by  Coley  et  al.  [1994]. 
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Our  magnetic  activity  analysis  study  showed  that  storm-time  events  affect  the 
neutral  winds  over  Arecibo  before  Arequipa.  This  analysis  also  indicated  that  the 
response  of  the  Arecibo  neutral  winds  to  disturbed  conditions  occur  within  3  hours  of 
the  storm  onset.  At  Arequipa,  we  showed  that  it  takes  longer  than  3  hours  after  the 
storm  onset  for  the  neutral  winds  to  respond  to  storm-time  effects.  We  require  a 
magnetic  activity  index  with  a  smaller  time  resolution  than  the  3-hour  Kp  to  provide  a 
more  precise  short-term  disturbed  analysis.  The  1-hour  AE  index  is  a  possible  candidate 
for  enhancing  the  short-term  disturbed  study.  The  main  drawback  is  the  limiting  of 
neutral  wind  statistics  because  thus  far  AE  index  has  routinely  been  calculated  only  up 


to  September  1987. 
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